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PREFACE 


The  High  Altitude  Effects  Simulation  (HAES)  Program  sponsored  by 
the  Defense  Nuclear  Agency  since  the  early  1970  time  period,  comprises 
several  groupings  of  separate,  but  interrelated  technical  activities,  e.g. 

ICECAP  (Infrared  Chemistry  Experiments  - Coordinated  Auroral  Program). 

Each  of  the  latter  have  the  common  objective  of  providing  information 
ascertained  as  essential  for  the  development  and  validation  of  pi^dictive 
computer  codes  designed  for  use  with  high  priority  DoD  radar,  communica- 
tions, and  optical,  defensive  systems. 

Since  the  inception  of  the  HAES  Program,  significant  achievement 
and  results  have  been  described  in  reports  published  by  DNA,  participating 
service  laboratories,  and  supportive  organizations.  In  order  to  provide  • 

greater  visibility  for  such  information  and  enhance  its  timely  applications, 
significant  reports  published  since  early  calendar  1974  shall  be  identified 
with  an  assigned  HAES  serial  number  and  the  appropriate  activity  acronym 
(e.g.  ICECAP)  as  part  of  the  report  title.  A complete  and  current  biblio- 
graphy of  all  HAES  reports  issued  prior  to  and  subsequent  to  HAES  Report 
Mo.  1,  dated  5 February  1974  entitled,  "Rocket  Launch  of  an  SWIR  Spectrometer 
into  an  Aurora  (ICECAP  72),"  AFCRL  Environmental'  Research  Paper  No.  466,  is 
maintained  and  available  on  request  from  DASIAC,  DoD  Nuclear  Information  and 
Analysis  Center,  816  State  Street,  Santa  Barbara,  California  93102,  Telephone 
(805)  965-0551. 

This  is  the  sixth  report  in  the  HAES  series  and  covers  technical  / 

work  performed  under  DNA  Contract  001-74-C-0143  from  1 January  1974  to 
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30  September  1974.  Po. tions  of  this  program  were  previously  reported  in 
DNA  Report  3297F  entitled:'  "Studies  of  Auroral  Simulation".  The  main 
purpose  of  the  work  is  to  understand  auroral  data,  especially  those 
obtained  under  the  ICECAP  program,  so  that  better  models  can  be  constructed 
to  describe  the  nuclear-disturbed  environment.  One  of  the  principal  tools 
used  for  this  purpose  is  the  ARCTIC  code  that  embodies  a wide  range  of 
phenomenon  including  excitation  of  atomic;.,  ionic,  and  molecular  states, 
chemical  reactions,  and  optical/IR  radiation  in  an  auroral  environment. 

The  present  report  includes  a description  of  some  calculations  pertinent 
to  auroral  events  and  comparisons  with  experimental  observations; 
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SECTION  1 
INTRODUCTION 


This  report  describes  the  continuing  work  on  auroral  simulation 
performed  from  1 January  1974  to  30  September  1974.  Related  work  on  the 
prior  contract  was  reported  in  Reference  1. 

The  main  purpose  of  these  studies  is  to  understand  auroral  data, 
especially  those  obtained  under  the  ICECAP  program,  so  that  better  models 
can  be  constructed  to  describe  the  nuclear-disturbed  environment. 

One  of  the  principal  tools  developed  for  use  in  these  studies  is 
the  ARCTIC  code,  described  in  Reference  1.  It  is  especially  suitable  for 
analyzing  data  from  multi-experiment  rocket  flights  that  provide  simul- 
taneous information  on  particle  fluxes,  chemical  species,  and  optical/ 
infrared  emissions  for  single  auroral  events.  Indeed,  it  is  in  this  capa- 
city that  the  code  has  been  largely  used  so  far. 

Section  2 below  contains  the  bulk  of  our  work  relating  to  use  of 
the  ARCTIC  code.  It  includes  (1)  the  detailed  energy  partition  by  primary 
electrons  incident  on  "normal"  air  as  predicted  by  the  current  version  of 
the  ARCTIC  code,  (2)  results  prepared  for  a code- comparison  meeting  at 
Lockheed,  Palo  Alto,  on  14  June  1974,  (3)  a detailed  comparison  with  ICECAP 
data  for  an  IBC  III  auroral  event  recorded  by  ground-  and  rocket-based 
instruments,  and  (4)  a comparison  between  calculated  and  observed  auroral 
emissions  from  02 ( 1 A). 
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Section  3 outlines  our  exploratory  efforts  to  account  for  certain 
long  wavelength  infrared  (LWIR)  features  observed  on  ICECAP.  It  also  gives 
some  preliminary  conclusions  with  respect  to  the  importance  of  aurorally- 
generated  extreme  ultraviolet  (EUV)  radiation  in  an  auroral  environment. 

Section  4 presents  the  results  of  our  study  of  atmospheric  heave 
under  intense  auroral  bombardment  conditions . 

Section  5 summarizes  our  conclusions  and  lists  some  specific 
recommendations. 

Three  appendices  are  included.  Appendix  A describes  a derivation 
of  the  spectral  distribution  of  chemiluminescent  emission  in  the  fundamental 
and  first  overtone  bands  of  NO  arising  from  the  reaction  N(ZD)  + 02  *►  NO  + 0. 
Appendix  B gives  corresponding  results  for  N0+  chemiluminescence  arising 
from  the  reactions  02+  + NO  ■+  02  + N0+  and  N2+  + 0 •*  N0+  + N.  Appendix 
C briefly  describes  the  structure  of  the  ARCTIC  code  and  includes  (on  micro- 
fiche) a complete  listing  of  the  program. 

Tentative  conclusions,  based  on  a study  of  the  "Murcray  phenom- 
enon", were  reported  separately  (Reference  2). 


SECTION  2 

APPLICATIONS  OF  THE  ARCTIC  CODE 


A detailed  description  of  the  ARCTIC  code  can  be  found  in  Reference 
1.  It  will  not  be  repeated  here,  although  a listing  of  the  code  and  * brief 
outline  of  its  structure  appear  in  Appendix  C.  Before  proceeding  with  the 
applications,  however,  the  following  remarks  are  appropriate. 

As  explained  in  Reference  1,  scattering  of  the  primary  electrons 
can  be  treated  using  two  different  versions  of  the  code.  The  first,  and 
simplest  version,  is  based  on  what  we  call  the  mean  scattering  approximation 
in  which  each  electron  is  assumed  to  have  the  same  range.  The  second  version 
called  scattering  with  range  variance,  recognizes  the  fact  that  there  is 
actually  a distribution  in  range  among  those  electrons  that  have  a common 
mean  range.  The  second,  and  more  accurate  version,  was  believed  to  be  con- 
siderably more  expensive  to  run  than  the  first  one.  We  have  since  found, 
however,  that  the  range -variance  version  is  much  faster  than  expected  (less 
than  a factor  of  2 more  than  the  simpler  version).  The  results  presented 
below  are  based  on  the  more  accurate  range-variance  version  of  the  code. 

This  section  presents  four  applications  of  the  code.  The  first, 
and  simplest,  is  a tabulation  of  species  populations  produced  in  "normal" 
air  by  the  total  stopping  of  primary  electrons.  Such  a tabulation  is  often 
needed  as  input  to  chemistry  and  optical/infrared  codes  used  for  weapon 
effects  studies  involving  the  deposition  of  x rays  or  beta  particles.  The 
second  is  a calculation,  for  the  case  of  a hypothetical  aurora,  performed 
for  comparison  purposes  at  a DNA- sponsored  meeting  at  Lockheed,  Palo  Alto, 
on  14  June  1974.  The  third,  and  major  application,  is  to  the  case  of  an 
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auroral  arc  measured  by  ground-  and  rocket-based  instruments  on  the  night 
of  March  27,  1975.  A detailed  comparison  with  the  data  available  to  us  is 
made,  and  conclusions  and  inferences  drawn  therefrom  are  given.  The  fourth 
application  is  to  auroral  enhancement  of  02  C1 A)  emission  at  1.27  pm 
compared  with  some  observations  reported  in  the  literature. 

SPECIES  POPULATIONS  IN  "NORMAL"  AIR  UNDER  ELECTRON  BOMBARDMENT 

The  partition  of  energy  into  states  of  excitation,  ionization,  and 
dissociation  was  determined  for  primary  electrons  stopped  in  a "normal"  air 
mixture  consisting  of  79%  N2  and  219«  O2 . Results  were  obtained  for 
primary  energies  between  0.17  eV  and  27  keV.  For  the  sake  of  compactness, 
however,  only  the  results  for  energies  of  0.92  keV  and  6 keV  are  presented 
here.  These  appear  in  Table  2-1. 

For  each  of  the  two  primary  enei-gies,  Table  2-1  lists  the  popu- 
lation of  the  species  produced  directly  by  the  primary  electron  as  well  as 
the  total  population,  including  contributions  from  all  generations  of 
secondary  electrons.  Also  listed  are  the  populations  produced  per  ion  pair. 

The  Rydberg  states  with  different  values  of  the  principal  quantum  number,  n, 

* 

are  lumped  together  for  each  of  the  series.  The  energies  of  the  individual 
states  fall  within  the  ranges  indicated  in  the  second  column  of  the  table. 

Calculated  values  of  the  average  energy  required  to  produce  an  ion 
pair  in  the  air  mixture,  as  a function  of  the  primary  electron  energy,  are 
shown  in  Figure  2-1.  They  are  in  excellent  agreement  with  limited  data 
obtained  at  energies  above  1 keV  with  betas  particles  in  air  (Reference  3) . 
Similar  calculations  have  been  performed  for  electrons  incident  in  pure  N2 
and  pure  02  and  the  results  (not  shown)  are  in  correspondingly  good  agree- 
ment with  the  limited  data. 


The  code  actually  computes  the  partition  into  states  with  different  n, 
but  for  the  sake  of  compactness  we  have  merged  the  results  in  Table  2-1 
for  each  of  the  Rydberg  series . 
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Table  2-1.  Species  populations  in  "normal" air  (79%  N2,  21%  02) 
produced  by  the  stopping  of  0.92  keV  and  6 keV 
electrons. 


Species 

State  Energy 
(eV) 

Primary 

Population 

Tot. 

Popula 

il 

tion 

Population 
Per  Ion  Pair 

0.92  keV 

6 keV 

0.92  keV 

6 keV 

0.92  keV 

6 keV 

n2(D 

a) 

0.29 

2.7 

2.7 

4.51 

2.92 

1.8 

1.7 

N2  (2 ) 

0.57 

9.4"1 

9.4"1 

1.6* 

l.O2  t 

6.3"1 

5.7”1 

N2(3) 

0.86 

3 . 8" 1 

3.8“ 1 

6.6 

4.21 

2.6-1 

2.4"1 

N2(4) 

1.13 

2.8"1 

2.8-1 

4.7 

3.01 

1.8"1 

1.7“* 

N2  (5) 

1.41 

1.8-1 

1.8"1 

3.1 

2.01 

1.2"1 

l.l-1 

N2(6) 

1.68 

1.2"1 

1.2-1 

2.1 

1.3* 

8.2"2 

7.4-2 

M7) 

1.95 

9.3"2 

9.3"2 

1.6 

l.O1 

6.3“2 

5.7“2 

n2 

(8) 

2.21 

5.5"2 

5.5“2 

9.4-1 

6.0 

3.7“2- 

3.4-2 

n2 

(A3ZJ) 

6.17 

5.2-1 

5.2-1 

5.1 

3.31 

2.0-1 

1.9-1 

n2 

(B3na) 

7.35 

6.6-1 

8.5"1 

4.5 

2.81  : 

1.8-J 

1.6-1 

n2 

(C3Hu) 

11.03 

2.6-1 

2.6“1 

1.7 

l.l1 

6.6-2 

6.3  2 

n2 

> 11.9 

3.1”2 

3.1”? 

2. 1-1 

1.3 

8. 2" 3 

7.4-3 

n2 

(a'n5) 

(2> 

8.55 

5.9”1 

1.0 

2.9 

1.81 

l.l-1 

l.O"1 

n2 

(3) 

12.25 

8.1"2 

1.6"1 

2.8-1 

1.7 

l.l-2 

9.7-3 

n2 

(b1^) 

(4) 

12.80 

1.9 

9.3 

2.8 

1.81  i 

l.l-1 

l.o-1 

n2 

(b-*T*) 

14.00 

8.7-1 

4.2 

1.2 

7.9 

4.7~2 

4.5"2 

n2 

[Ry(x] 

P 

13.01  14.84 

2.6 

1.31 

3.6 

2.41 

1.4-1 

1.4"1 

n2 

:Ry(A)?» 

13.19  -»•  15.86 

7.5-1 

3.6 

1.0 

6.3 

3.9-2 

3.9"2 

n2 

15.75  -*•  17.95 

8.1-2 

4.2"1 

l.O"1 

6.8-1 

3.9-3 

3.9"3 

n2 

:Ry(D)] 

15.71  - 20.87 

1.4"1 

7.2"1 

1.8-1 

1.2 

7.0-3 

6.9"3 

n2 

•Ry(c)] 

20.79  - 22.83 

3.4-2 

1.9-1 

4.0-2 

2.7"1 

1 .6“ 3 

1.5"3 

N(4S) 

(7) 

- 

3.7 

1.81 

5.8 

3.81 

2.3-1 

2.2"1 

N(2D) 

- 

3.7 

1.91 

5.4 

3.71 

2.1-1 

2. 1"1 

N(2P) 

1.1 

6.7 

1.6 

1.21 

6.3~2 

6.7-2 

Nt 

(x2ng) 

(8) 

15.58 

7.2 

3.61 

9.0 

5.91 

3.5”1 

3.4" 1 ! 

Nt 

(A2nSl 

16.73 

3.6 

1.91 

4.4 

2.91 

1.7-1 

1.7"1 

Nt  B2tf 

18.75 

1.6 

8.5 

1.8 

i.31 

7.0-2 

7.4"2 

Nt 

(D*n5-! 

22.00 

7.1 -1 

4.1 

8. 1"1 

5.6 

3.2"2 

3.2-2 

Nt(C2ESI 

23.60 

7.3”1 

4.3 

8. 1"1 

5.8 

3.2-2 

3.3"2 

N* 

(3P) 

0) 

• 

1.9 

1.21 

2.1 

Wm 

£ MB| 

N+(2D) 

• 

7.9-1 

5.7 

8.8"1 

■SSI 

mssm 

N+(JS) 

- 

2.4- 1 

1.8 

2.7"1 

am 

turn 

1 
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Table  2-1  (Continuer) . Species  populations  in  "normaV'air  (79%  N2.  21%  02) 

produced  by  the  stopping  of  0.92  key  and  6 keV 
electrons. 


Species 

State  Energy 
(eV) 

Primary 

Population 

Total  j 

Population  j 

Population 
Per  Ion  Pair 

0.92  keV 

6 keV 

0.92  keV 

6 keV 

0.92  keV 

0 2 ( l/  ) .... 

-0.25 

5.6 

5.6 

1.1 2 

7.22 

4.3 

4.1 

9.9  eV  peak*10 

9.90 

6.7-2 

2 . 4“ 1 

1 . 2“ 1 

7 . 5" 1 

4.7"3 

4. 3" 3 

(MaMg) 

0.98 

6.4-1 

6.4"1 

9.9 

6 . 4 1 

3.9"1 

3.7"1 

Mb'IS)  „„ 

1.64 

1.5-1 

1.5"1 

2.2 

1 . 4 1 

8.6"2 

8.0"2 

°2(A* 

4.50 

6.2-1 

7.5-1 

7.3 

4.71 

2.9"1 

2.7"1 

02(b3eg)  U2) 

6.10 

1.0 

3.8 

2.4 

1.51 

9.4"2 

8.6"2 

OaCRytX)],  , 

8.41  ->  11.22 

3.5-1 

1.5 

6.5"1 

4.1 

2.5"2 

2.3“2 

02[Ry{a)]'1!l 

12.70  - 15.25 

1.7-1 

8.4'1 

2.4“ 1 

1.6 

9.4"3 

9. 1"3 

02[Ry(A)] 

13.25  -v  16.02 

1.5-1 

7 . 3“ 1 

2.0"1 

1.3 

7.8"3 

7. 4" 3 

02[Ry(bj] 

15.65  + 17.47 

6. 4- 2 

3.3"1 

8.2"2 

5.4"1 

3.2“3 

3.T3 

02[Ry(B»c)](w) 

19.60  -v  22.15 

2.5"2 

1.4"1 

2.9"2 

2.0"1 

1J"3 

l.l"3 

0(3P) 

- 

3.9 

1.41 

1-.7-1 

i.r- 

6.6"1 

6.3"1 

0(‘D) 

- 

1.0 

3.8 

2.4 

1.51 

9.4"2 

8.6"2 

12.10 

7.3"1 

3.6 

9.9"1 

6.5 

3.9“2 

3.7"2 

ot(a4nu) 

16.10 

7.0"* 

3.6 

8.5"1 

5.6 

3.3"2 

3.2“2 

o*(A2nu) 

16.90 

6.8"1 

3.6 

8.r! 

5.4 

3.2“2 

3. 1“2 

0t(b4Zg) 

18.20 

5.6"1 

3.0- 

6.6"1 

4.4 

2.6"2 

2.5“2 

0t(B,cJ  <15> 

23.00 

2.4-1 

1.4 

2.6"1 

1.9 

1.0"2 

l.l"2 

o^rs) 

- 

1.1 

5.4 

1.2 

7.7 

4.7"2 

4. 4"  2 

0+(2D) 

6.2"1 

3.2 

6 . 7” 1 

4.4 

2.6"2 

2.5“2 

(1) 

(2) 
(3) 

M 

(5) 

(6) 


(?) 

(8) 

(9) 

110) 

(111) 

(12) 


N2(v)  refers  to  the  v*n  vibrational  state  of  (X 1 S) . 

18%  predissociation  into  N('*S)  + N(4S)  assumed. 

The  lowest  Rydberg  state  of  N2. 

69%  predissociation  into  N("S)  + N(2D)  assumed. 

Ry(X)  means  the  Rydberg  states  (n  = 3,4,5+---)  converging  on  N2+(X2Zq); 
similarly  Ry (A)  means  the  Rydberg  states  converging  on  N2+  (A2nu),  etc. 
For  the  Rydberg  states  Ry[A(n=5),B,D,C],  the  numbers  in  the  table  are 
25%  of  the  total  excitations  to  the  respective  states-  The  remaining 
75%  are  assumed  to  be  autoionized  and  have  been  added  to  N2+(X'jig). 

The  nitrogen  atoms  arise  from  predissociation  and  from  disscc  ative 
excitation  and  dissociative  ionization. 

Includes  contributions  from  autoionization  of  Rydberg  stares.  See  note 
(6)  above). 

The  atomic  ions  arise  from  dissociative  ionization. 

Unassigned. 

Assumed  to  be  100%  dissociated  into  0(3P)  +0(3P). 

Assumed  to  be  100%  dissociated  into  0(3P) -*  0( *D). 
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Table  2-1  (Continued).,  Specie?  popuMtions  in  "normal" air  (79%  N2»  21%  02) 

produced  by  the  stopping  of  0.92  keV  and  6 keV 
electrons. 


(13>  For  the  Rydberg  states  Ry(a, — ),  the  numbers  in  the  table  are  50% 
of  the  total  excitations  to  the  respective  states.  The  remaining 
50%  are  assumed  to  be  autoioni2ed  and  have  been  added  to  O2 (X2IIg) . 
(14>  Ry(B,c)  is  a blend  of  the  Rydberg  states  converging  on  the  B2E<j 
and  c4ZQ  states  of  02+. 

(15)  A blend  of  the  states  B2Zg  and  c4Z[j. 


It  should  be  pointed  out,  however,  that  the  agreement  with  data 
is  partly  forced  because  the  unknown  autoionization  factors  have  been 
adjusted  to  improve  the  agreement.  Furthermore,  a large  number  of  the 
excitation  cross  sections  are  not  well  known  and,  in  some  cases,  especially 
for  the  Rydberg  states,  are  little  more  than  guesses.  Changes  in  the  rela- 
tively large  cross  sections  for  the  Rydberg  states  can  have  an  appreciable 
effect  c the  calculated  values  for  "eV  per  ion  pair". 

TEST  CASE  '( HYPOTHETICAL  AURORA) 

For  the  code- comparison  meeting  at  Lockheed,  two  test  problems 
were  run  using  the  ARCTIC  code.  The  first  one,  described  here,  involved  a 
hypothetical  electron  flux  of  exponential  type.  The  second  one,  described 
in  detail  in  the  next  subsection,  was  based  on  a flux  spectrum  measured 
under  the  ICECAP  program. 

The  first  problem  was  essentially  as  follows.  Given  an  incident 
spectral  flux  of  primary  electrons: 

flux  = e”^10  (electrons  cm"2  sec"1  ster"1  keV”1) 
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ted  with  the  ARCTIC  code. 


at  an  altitude  of  250  km,  isotropic  downward  with  a vertical  magnetic  field, 
determine  the  energy  deposition  rate  as  a function  of  altitude,  and  the 
volume  production  rate  of  selected  states.  The  atmosphere  specified  was 
CIRA  1965,  model  5,  hour  0. 

Figure  2-2  shows  the  calculated  energy  deposition  rate  profile. 
Shown  for  comparison  (dashed  curve)  are  the  results  obtained  with  the 
Lockheed  (LMSC)  code.  They  are  seen  to  be  in  substantial  agreement  vith 
those  from  the  ARCTIC  code.  A study  of  the  differences  has  shown  tl  ac  they 
arise  almost  entirely  from  differences  in  the  assumed  elastic  scattering 
cross  sections.  Results  obtained  with  the  NRL  and  Visidyne  codes  were  at 
variance  with  each  other  as  well  as  with  those  in  igure  2-2.  They  are  not 
shown  here. 

Production  rates  for  about  200  different  excited  states  of  N2, 

02,  N,  and  0 were  obtained.  Figures  2-3  to  2-6  show  a few  selected 
results  from  the  list.  Comparison  of  these  results  with  those  from  other 
codes  was  not  made, 

AURORAL  ARC  OF  27  MARCH  1973  (BLACK  BRANT  18.205-1) 

Preliminary  Discussion 

A considerable  fraction  of  our  effort  has  been  devoted  to  a study 
of  the  aurorr.l  event  of  27  March  1973,  observed  in  the  vicinity  of  Poker 
Flat,  Alaska.  Details  of  the  fvent  (a».  arc  with  brightness  coefficient  at 
least  11+)  were  recorded  by  instruments  on  board  a Black  Brant  rocket  and 
by  various  ground-based  instruments  including  those  located  at  Ester  Dome 
(EDO)  and  Fort  Yukon  (FYU):,  Alaska.  All  of  the  data  obtained  with  these 
instruments  are  not  yet  reduced  and  available.  The  analyses  and  results 
reported  below  are  based  on  the  limited  amount  of  preliminary  data  that 
was  made  available  to  us.  We  believe  that  a better  model  to  describe  the 
event,  especially  the  auroral  environment  outside  the  main  arc,  can  I z 
constructed  when  all  of  the  evidence  is  at  hand. 
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Incident  flux  (250  km)  = e-^' ,0( electrons  cm-2secT1ste 


description  of  incident  flux,  see  Figure  Z-3 


description  of  incident  flux,  see  Figure  2-3 


For  description  of  incident 
see  Figure  2-3. 


In  this  subsection  we  will  first  give  a brief  overview  of  the 
auroral  event,  together  with  a description  of  the  types  of  data  that  were 
taken  on  it.  Next  we  proceed  to  set  the  stage  for  actual  calculations  with 
the  ARCTIC  code  by  describing  the  environment  into  which  the  electron  flux 
in  the  main  ai'c  is  assumed  to  be  deposited.  This  description  includes  a 
consideration  of  the  relatively  soft  spectral  flux  of  electrons,  that  we 
will  call  a "drizzle",  tha“  is  believed  to  have  bombarded  the  air  over  a 
large  region  for  some  time  prior  to  the  rocket  launch.  We  next  go  on  to 
deduce  some  geometrical  and  dynamical  characteristics  of  the  arc  itself, 
and  present  the  data  used  to  describe  the  incident  spectral  flux.  The  cal- 
culations are  then  described  and  the  results,  together  with  data  comparisons, 
are  presented.  Finally,  some  conclusions,  based  on  the  comparisons,  are 
drawn. 

Qualitative  Description  of  Event 

A Black  Brant  rocket  (18.205-1)  was  launched  from  Poker  Flat  (EF) , 
Alaska  on  27  March  1973  at  0957:45  UT  into  a bright  auroral  arc.  The 
geometrical  arrangement,  including  the  projection  of  the  rocket  trajectory 
in  the  magnetic  meridian  plane  and  the  approximate  position  of  the  arc  at 
91  sec  after  rocket  launch,  is  shown  in  Figure  2-7.  Details  concerning  the 
rocket  attitude  and  trajectory  and  the  on-board  instrumentation  can  be 
found  in  Reference  4.  The  geometry  depicted  in  Figure  2-7  is  based  on 
information  supplied  to  us  by  Photometries  (Reference  5).  For  altitudes 
above  about  70  km,  the  rocket  axis  was  nearly  vertical,  being  inclined 
3 degrees  south  (magnetic)  from  the  zenith,  and  10  degrees  from  the 
magnetic  field  lines. 

Although  the  rocket  was  programmed  to  enter  the  arc  at  100-km 
altitude,  it  actually  entered  between  105  and  110  km  and  overflew  the  region 
of  maximum  energy  deposition  (Reference  6).  the  on-board  instrumentation 
that  is  relevant  to  this  report  included:  (1)  an  electrostatic  analyzer  to 
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measure  the  differential  electron  spectra  over  the  energy  range  3 to  30  keV, 
(2)  a scintillator  to  measure  the  energy  deposited  by  electrons  with  energies 
greater  than  about  3 keV,  (3)  a positive  ion  and  neutral  mass  spectrometer, 

(4)  plasma  frequency  and  Langmuir  probes  to  measure  the  electron  density, 

(5)  a vertical-viewing  photometer  at  3914  X (in  the  rocket  nose),  (6)  a 

vertical- viewing  dual  channel  radiometer  at  4.3  and  5.4  ym,  (7)  a vertical- 
viewing CVF  spectrometer  (8)  a side- looking  (80°  from  vertical)  dual  channel 
radiometer  at  2.7  and  5.4  ym,  and  (9)  five  side- looking  (80°  from  vertical) 

0 0 0 0 o 

optical  photometers  at  3466  A,  3800  A,  3914  A,  5199  A,  and  5577  A. 

O O 

Ground-based  instrumentation  included  4278  A-  and  5577  A-scanning 
photometers  (scanning  in  the  magnetic  meridian  plane)  at  Ester  Dome  and 
Fort  Yukon  (see  Figure  2-7)  as  well  as  all-sky  cameras  at  both  of  these 
locations. 


The  general  position  and  motion  of  the  arc,  as  deduced  from  the 
scanning  photometer  data  at  Ester  Dome  and  Fort  Yukon  (Reference  5)  in  the 
time  span  from  about  10  minutes  before  to  7 minutes  after  rocket  launch  is 
as  follows.  At  -10  minutes  the  arc  was  nearly  overhead  at  Ester  Dome. 

It  subsequently  moved,  rather  erratically,  northward  until  it  reached  Poker 
Flat  at  about  -9  minutes.  It.  remained  slightly  north  of  that  position 
until  about  -6.5  minutes  when  it  moved  northward  again  and  reached  its 
most  northerly  position  (about  halfway  between  Ester  Dome  and  Fort  Yukon) 
at  about  -5  minutes.  It  then  proceeded  to  move  south  until  about  95  sec 
after  rocket  launch  (approximately  the  time  the  rocket  entered  the  main 
arc  — when  the  arc  was  ~ 48  km  north  of  Poker  Flat)  when  it  moved  northward 
again.  An  erratic  movement  northward  continued;  until  about  +4.5  minutes, 
at  which  time  the  arc  was  70  km  north  of  Poker  Flat,  and  a generally  south- 
ward movement  then  ensued. 

The  foregoing  description  should  serve  to  indicate  that  a large 
air  volume,  including  the  measurement  space  itself,  was  bombarded  sporadically 
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by  the  arc  flux  for  some  time  prior  to  the  rocket  measurements.  In  fact, 
the  point  where  the  rocket  first  penetrated  the  arc  had  been  irradiated 
previously  at  about  -5.5  minutes  as  the  arc  moved  northward  through  it. 

Exo-Arc  Environment 

It  is  clear  from  the  above  discussion  that  the  photometer-  and 
radiometer-sight  paths  traversed  regions  outside  the  arc  that  had  previously 
been  irradiated  in  a complex  and  ill-defined  manner.  In  order  to  make 
meaningful  comparisons  between  ARCTIC-code  calculations  and  the  data,  we  have 
attempted  to  construct  a mean  exo-arc  environment  that  may  suffice  for  the 
purpose.  This  model  environment  will  now  be  described. 

'Hie  assumed  initial  neutral  atmosphere  (prior  to  any  irradiation) 
is  shown  in  Table  2-2. 


Table  2-2.  Assumed  initial  nighttime  neutral  atmosphere. 


Quantity 

Altitude  Range 
(km) 

Values 

[N2],  [Oil,  [0],  T 

J 80  < hi  < 120  * 
(l20  < h < 800 

1965  CIRA  mean,  p . 1 2 
1965  CIRA  mean,  mod. 5,  hr.0 

[03] 

> 30 

4.5  x 1012  exp  [ (30-h )/5 . 42] 

[C02] 

< 100 
> 100 

3 x 10"4  [molecules] 

3 x no" 4 [molecules]exp[(100-h)/17.4] 

< 80 

4 x ]o" 8 [molecules] 

[NO] 

80  < h < no 

4 x ]0"8  [molecu1es]exp[(h-80)/4.29] 

> no 

4 x TO"8  [molecules]exp[6.99+ 
(110-h)/43.5] 

* h is  altitude  in 

kilometers. 
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The  vibrational  states  of  N2  are  assumed  to  be  Boltzmann  distri- 
buted initially,  and  the  electron  temperature  is  equated  to  that  of  the 
neutral  gas. 

An  error  in  the  assumed  vertical  distribution  of  the  major  species 
can  lead  to  several  sources  of  error  in  the  calculated  chemical  and  optical 
environment.  Included  among  these,  for  example,  are  errors  in  the  altitude 
profile  of  the  energy  deposition  rate,  in  the  energy  partition  as  a function 
of  altitude,  and  in  the  rate  at  which  chemical  reactions  can  proceed.  Un- 
fortunately, neutral  mass  spectrometer  data  for  this  event  are  not  currently 
available  to  us  and,  consequently,  the  accuracy  of  our  assumed  profiles  for 
the  major  species  is  certainly  in  question. 

Considerable  evidence  is  available  to  suggest  that  a weak  particle 
bombardment  is  present  in  the  auroral  oval  at  all  times  (References  7,  8,  9). 
In  particular.  Walker  et  al.,  (Reference  8)  have  fit  one  set  of  observations 
with  a bombarding  flux  of  the  form  E exp(-E/0.6),  isotropic  over  the  lower 
hemisphere.  Using  these  results,  Jones  and  Rees  (Reference  10)  have  calcu- 
lated an  "ambient"  high- latitude  nocturnal  ionosphere.  We  have  used  their 
results  to  define  relative  concentrations  for  N0+,  02+,  and  0+,  and  have 
assumed  an  "ambient"  ionospheric  model  by  normalizing  these  results  to 
those  of  Walker  (Reference  11).  Concentrations  of  all  other  species  are 
initially  set  to  zero. 

The  foregoing  remarks  describe  the  "ambient"  high-latitude  atmo- 
spheric model  assumed.  Into  this  environment  must  now  be  superimposed  the 
modifications  that  arise  from  the  intense,  but  localized,  electron  flux  into 
the  auroral  arc,  as  well  as  that  from  the  associated  drizzle  precipitation 
that  extends  well  beyond  the  arc. 

The  work  of  Frank  and  Ackerson  (Reference  9)  shows  that  the 
particle  flux  outside  an  arc  is  associated  in  a complex  way  with  that  of 
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the  arc  itself.  It  is  generally  not  a constant  in  space  but  shows  con- 
siderable structure  of  the  "inverted-V"  type.  That  is,  the  hardness  and 
strength  of  the  particle  spectrum  tends  to  fall  off  with  distance  north 
and  south  from  the  arc  proper.  The  arc  and  its  surrounding  flux  (that 
may  extend  for  hundreds  of  kilometers  north  and  south  of  the  arc)  evidently 
constitute  a unified  event.  For  convenience  we  shall  refer  to  the  localized 
region  of  intense  precipitation  as  the  "main  arc",  and  the  larger  region  of 
complex  but  weaker  precipitation  as  the  "drizzle"  region. 

Since  the  main  arc  generally  moves  within  the  drizzle  region  we 
shall,  for  calculational  purposes,  assume  that  the  two  regions  are  essen- 
tially independent  and  that  the  drizzle  flux  serves  to  determine  the  back- 
ground environment  into  which  the  main  arc  flux  is  deposited. 

The  only  particle  flux  data  available  for  the  drizzle  region  are 
those  obtained  over  that  part  of  the  rocket  trajectory  outside  the  main 
arc  (see  Figure  2-7).  Samples  of  such  data  (Reference  6),  corresponding 
to  various  times  after  rocket  launch,  are  shown  in  Figure  2-8*  For  calcu- 
lational simplicity,  we  have  adopted  an  "average"  drizzle  flux  spectrum 
shown  by  the  dashed  curve  in  Figure  2-8,  and  assume  that  the  precipitation 
is  isotropic  (in  pitch  angle)  in  the  downward  hemisphere.  The  adopted 
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spectrum  was  normalized  so  as  to  provide  agreement  with  the  measured  4278  A 
zenith  brightness  at  Fort  Yukon.  This  flux  is  assumed  to  be  constant  in 
space  and  time  preceding  and  during  the  rocket  flight,  although  the 
assumption  is  clearly  at  variance  with  the  general  findings  of  Frank  and 
Ackerson  mentioned  above.  For  this  reason  we  cannot  hope  to  achieve  precise 
agreement  with  the  data  from  some  of  the  side- looking  instruments  on  board 
the  rocket  (especially  the  5199  A and  3466  A emissions  from  NO  D)  and 
N(2P),  respectively)  that  depend  rather  sensitively  on  details  of  the  drizzle. 

★ 

These  spectra  show  many  of  the  characteristics  of  an  intense  "inverted  V" 
structure,  particularly  the  progressive  softening  and  weakening  of  the 
spectrum  as  a function  of  distance  from  the  main  arc. 
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Figure  2-8.  Assumed  drizzle  spectrum  compared  with  rocket  (Black  Brant 
18.205-1)  data  outside  arc. 
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Rather,  in  these  cases,  the  object  is  to  account  for  the  general  magnitude 
of  the  observed  emissions  which,  in  the  absence  of  an  assumed  drizzle,  we 
would  not  otherwise  be  able  to  do. 

In  addition  to  providing  a drizzle  spectrum,  it  is  also  necessary 
to  specify  the  time  over  which  it  has  operated  and  its  spatial  extent. 

Little  information  is  available  on  the  effective  duration  of  the  drizzle 
bombardment,  although  data  from  the  side- looking  photometers  on  board  the 
rocket  provide  some  guide  as  to  its  spatial  extent.  With  respect  to  the 
bombardment  time,  however,  we  note  that  for  some  ten  days  preceding  the 
event,  significant  magnetic  disturbances  with  indices  ranging  from 

4 to  7 occurred  (Reference  12).  In  addition,  ground-based  photometer  data 
(Reference  5)  show  that  a drizzle  bombardment  probably  occurred  for  at 
least  Vi  hours  before  rocket  launch. 

In  order  to  determine  the  effect  on  species  concentrations  due  to 
a long  term  drizzle  precipitation,  we  ran  the  ARCTIC  code  with  the  drizzle 
spectrum  shown  in  Figure  2-8  (dashed  curve)  for  times  to  106  sec.  Figures 
2-9  and  2-10  show  the  resulting  profiles  of  0?,  0,  and  NO  for  selected 
drizzle  duration  times.  The  results  clearly  show  an  NO  buildup,  especially 
for  106  sec  (11.5  days),  along  with  an  increase  in  [0]  and  a decrease  in 
[02].  Since  the  effects  due  to  diffusion,  transport,  and  daytime  solar 
radiation  are  not  included  in  the  calculation,  the  late-time  results  are 
not  realistic.  However,  they  do  serve  to  indicate  the  type  of  chemical 
effects  that  can  arise  from  prolonged  bombardment. 

Large  enhancement  of  the  NO  concentration  in  an  auroral  environ- 
ment has  been  reported  by  Zipf  (Reference  13),  and  has  also  been  inferred 
from  numerous  ion  mass  spectrometer  measurements  (eg.,  References  14,  15). 

In  fact,  for  the  auroral  event  considered  here,  we  have  inferred  an  enhanced 
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NO  profile  using  the  positive  ion  mass  spectrometer  data  of  Narcisi  (Refer- 
ence 16)  and  a procedure  similar  to  that  described  by  Swider  and  Narcisi 
(Reference  15).  Because  of  rocket  spin  modulation,  the  mass  spectrometer 
data  show  large  fluctuations  (see  Figure  2-39),  and  the  inferred  NO  profile 

reflects  this  ambiguity.  However,  Figure  2-11  shows  an  NO  profile  that  is 
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not  inconsistent  with  the  mass  spectrometer  data.  Shown  for  comparison  is 
a measured  low  latitude  profile  (References  17,  18).  The  comparison 
suggests  that  the  NO  concentration  in  the  vicinity  of  the  auroral  arc 
was  greatly  enhanced  above  its  low  latitude  values,  at  least  below  about 
120  km,  and  attains  values  ~ 2 x 109  cm"3  between  100-  and  105-km  altitude. 

In  preparing  an  environment  in  which  to  deposit  the  main  arc  flux, 
we  have,  somewhat  arbitrarily,  assumed:  (1)  the  drizzle  flux  is  turned  on 
30  minutes  prior  to  rocket  launch,  (2)  all  concentrations,  with  the  exception 
of  NO,  are  initially  specified  by  the  high-latitude  "ambient"  model 
described  earlier,  (3)  the  NO  profile  is  given  initially  by  the  enhanced 
values  shown  in  Figure  2-11. 

Unfortunately,  as  indicated  earlier,  errors  in  the  assumed  con- 
centrations of  02,  0,  and  NO  can  lead  to  significant  differences  between 

the  computed  and  measured  environment.  For  example,  the  ion  concentration 
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ratios,  especially  NO  /O2  , are  sensitive  to  the  NO  profile,  and  detailed 
agreement  between  calculations  and  observations  should  not  be  expected  in 
the  absence  of  good  data  on  [NO] . 

With  the  foregoing  assumptions,  the  perturbed  atmosphere  (result- 
ing from  30  minutes  of  drizzle  bombardment):  into  which  the  main  arc  flux 
is  deposited,  is  shown  in  Table  2-3. 

Actuallv,  the  inferred  values  for  NO  are  higher  above  120  km  than  are 
shown  xn  Figure  2-il.  For  simplicity,  we  have  extrapolated  the  enhanced 
profile  above  120  km  until  it  intersects  the  low  latitude  curve. 
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CONCENTRATION  (crrT3) 


Table  2-3.  Perturbed  atmosphere  following  30  minutes  of  "drizzle"  bombardment. 
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Table  2-3  (Continued).  Perturbed  atmosphere  following  30  minutes  of  "drizzle"  bombardment. 
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As  to  the  spatial  extent  of  the  drizzle  region,  we  have  made 
estimates  based  on  comparisons  between  calculations  and  data  from  the  side- 
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looking  rocket-borne  photometers  and  the  ground-based  4278  A instruments. 

The  result  of  these  comparisons,  as  shown  in  Figure  2-7,  suggests  a drizzle 
region  extending  about  310  km  south  and  540  km  north  from  the  main  arc,  for 

a total  north-south  extent  of  about  850  km.  The  reasonableness  of  this 

conclusion  may  be  judged  by  comparing  it  with  the  previously  cited  work  of 
Frank  and  Ackerson  (Reference  9)  on  other  auroral  events.  They  show  three 
sets  of  data  (Plates  4a,  4b,  and  5a)  corresponding  to  indices  of  4-, 

4,  and  4+.  For  these  events,  the  north-south  extent  of  the  associated 
drizzle  was  800,  1200,  and  1400  km,  respectively.  At  the  time  of  the  rocket 
flight  considered  here  (27  March  1973) , a value  of  about  4 appears  to 

be  appropriate  (Reference  12).  The  inferred  extent  of  our  drizzle,  therefore, 
is  not  inconsistent  with  these  other  auroral  data. 

Arc  Geometry  and  Dynamics 

To  proceed  farther,  the  geometry  and  velocity  of  the  main  arc 

must  be  determined  as  well  as  the  relationship  between  the  arc  and  the 

rocket-  and  ground-based  instruments. 

Visual  inspection  of  the  flight  records  of  the  electrostatic 
analyzer  reveals  several  important  ;points  concerning  the  geometry.  At 
about  91  sec  the  rocket  entered  a region  of  intense  particle  precipitation 
and  then  left  it  at  about  121  sec.  Two  typical  plots  of  the  measured  flux 
spectrum  (Reference  6)  during  this  time  span  are  shown  in  Figure  2-12. 

The  time  required  to  enter  and  leave  this  region  was  very  short,  indicating 
sharp  spatial  boundaries. 

Knowing  the  altitude  at  which  the  rocket  entered  and  left  the 
arc,  we  can  then  determine  the  horizontal  width,  AW,  of  the  arc  (north-south 
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dimension)  if  the  horizontal  component  of  the  relative  velocity  between 
the  arc  and  the  rocket  in  the  magnetic  meridian  plane  is  known.  In 
particular, 

121 

AW  = Ah/tan  6 + J Vr  dt  (2-1) 

91 

where  Ah  is  the  altitude  difference  between  the  points  of  rocket  entry 
to  and  exit  from  the  main  arc,  6 is  the  magnetic  dip  angle  (77°),  and 
is  the  horizontal  component  of  the  relative  velocity  in  the  magnetic 
meridian  plane.  From  the  rocket  trajectory  data  supplied  to  us  (Reference 
6),  Equation  2-1  becomes,  for  constant  V^, 

AW  a 7.1  + SO.  V (km)  . (2-2) 

An  estimate  of  the  relative  velocity,  V , can  be  obtained  using 
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radiance  data  from  the  upward-viewing  3914  A photometer  in  the  rocket  nose. 

This  is  shown  in  Figure  2-13  along  with  several  simulations  of  this  emission 

for  rocket  altitudes  below  110  km.  Inspection  of  Figure  2-13  shows  that  the 

slope  of  the  calculated  radiance  values  below  110  km  is  very  sensitive  to 

the  assumed  relative  velocity..  We  have  concluded  that  a value  V *s  - 0.07 
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km  sec  1 best  fits  the  data.  This  implies  an  average  arc  velocity  north- 
ward of  0.61  km  sec"1,  slightly  faster  than  that  of  the  rocket  itself  and, 
from  Equation  2-2,  a horizontal  arc  width  of  5 km. 

It  should  be  emphasized  that  the  deduced  arc  motion  is  valid, 
at  best,  only  over  a very  limited  time  span.  In  fact,  the  general  descrip- 
tion given  earlier  of  the  arc  motion  prior  to- rocket  launch  indicated  that 
the  main  arc  was  actually  moving  south  just  prior  to  the  time  interval  con- 
sidered here.  However,  for  the  time  period  considered,  the  deduced  arc 
velocity  is  in  approximate  agreement  with  the  ground-based  photometer 
data  (Reference  5) . 
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Figure  2-1 


To  summarize,  the  model  of  the  auroral  arc  system  we  have  devised 
is  shown  in  Figure  2-7.  It  consists  of  a main  arc,  approximately  5 km 
wide  (north-south) , moving  0.61  km  sec-1  to  the  north,  imbedded  in  a drizzle 
region  that  extends  about  850  km  in  the  magnetic  meridian  plane.  The 
absolute  position  of  the  main  arc  is  determined  from  a knowledge  of  the 
rocket  location  when  it  entered  the  main  arc.  For  the  southern  magnetic 
boundary  of  the  arc  at  108. 7- km  altitude,  this  corresponds  to  a displace- 
ment from  Ester  Dome  of  approximately  80  km.  Although  the  drizzle  is 
probably  of  the  "inverted-V"  type,  we  have  assumed  it  to  be  constant  over 
the  entire  region. 

Electron  Flux  in  Main  Arc  and  Energy  Deposition 

Consider  now  the  electron  flux  into  the  main  arc.  We  have 
already  seen  (Figure  2-12)  two  measured  particle  spectra  for  this  event. 
Visual  inspection  of  the  flight  records  at  the  Air  Force  Cambridge  Research 
Laboratories  (AFCRL)  reveals  that  the  spectra  were  nearly  constant  while 
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the  rocket  remained  in  the  arc.  However,  the  ground-based  4278  A and  5577 
data  show  a slow  decrease  in  intensity  between  +100  sec  after  rocket  launch 
and  termination  of  the  flight  at  +413  sec  (Reference  51),  thus  suggesting 
at  least  a weak  time  variation  of  the  electron  source.  We  shall  ignore 
any  such  variations. 

The  simplest  assumption  concerning  the  pitch-angle  distribution 
is  that  it  is  isotropic  over  the  downward  hemisphere.  The  limited  data  in 
this  respect  (Reference  6)  indicate  that,  at  least  over  the  pitch-angle 
range  from  35°  to  55°,  this  assumption  is  valid.  However,  the  magnitude 
of  the  spectral  flux  shown  in  Figure  2-12,  based  on  this  assumption,  leads 
to  discrepancies  with  the  optical  data.  In  fact,  if  the  magnitude  of  the 
particle  data  is  reduced  by  a factor  of  2 to  3,  better  agreement  with  the 
upward-looking  rocket-based;  3914  A and  ground-based  42=78  X radiances  appears 
to  be  achieved.  Other  possible  alternatives  are  as  follows.  One  is  to 
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assume  an  arc  width  of  only  2 km.  This  would  then  imply  (Equation  2-2), 
however,  an  arc  velocity  of  0.71  km  sec’1  that  is  clearly  at  variance  with 
the  data  shown  in  Figure  2-13.  Other  possibilities  are  that  the  pitch- 
angle  distribution  is  non  isotropic  or  that  the  photometer  data  are  in 
error.  The  recent  work  of  Evans  (Reference  19),  in  fact,  shows  that 
precipitating  electrons  usually  exhibit  a complex  pitch-angle  distribution 
ranging  from  one  that  is  primarily  field  aligned  to  one  that  is  isotropic. 

In  order  to  proceed  with  an  ARCTIC  code  calculation,  a definite 
decision  on  what  to  input  for  the  particle  flux  had  to  be  made.  The  weight 
of  evidence  at  the  time  seemed  to  suggest  that  wo  should  use  the  particle 

spectra  shown  in  Figure  2-12,  reduced  by  a factor  of  2.  This,  indeed,  is 
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what  we  did.  In  so  doing,  greater  weight  was  thus  given  to  the  3914  A 
rocket  data  shown  in  Figure  2-13.  This  was  done  because  the  ground-based 
data  seemed  to  be  in  support  of  that  from  the  rocket-based  photometer. 

The  assumed  arc  flux,  incident  at  an  altitude  of  200  km,  is  shown 
in  Figure  2-14.  It  is  one  half  of  the  115  sec  flux  shown  in  Figure  2-12. 

No  data  at  energies  below  3 keV  are  available,  so  we  have  arbitrarily  ex- 
trapolated the  spectrum  at  3 keV  down  to  zero  energy,  although  this  proce- 
dure is  hazardous.  It  should  be  pointed  out,  in  fact,  that  most  of  the 
energy  deposited  above  about  125  km  is  by  electrons  with  energy  < 3 keV.. 
Consequently,  in  the  absence  of  measured  particle  data  for  this  energy 
region,  results  based  on  the  assumed  extrapolated  spectrum  are  definitely 
suspect.  Indeed,  as  we  will  see  later,  comparisons  between  our  calculated 
results  and  the  data  suggest  that  the  assumed  spectrum  (Figure  2-14)  below 
3 keV  is  probably  an  underestimate,  while  that  shown  in  Figure  2-8  for  the 
drizzle  region,  is  probably  an  overestimate.  Indirect  evidence  of  an  en- 
hanced arc  flux  below  3 keV  can  be  seen  from  the  data  curve  in  Figure  2-13 
(also  Figure  2-23)  where  a change  of  slope  occurs  at  about  125  km. 

All  results  shown  below  in  this  subsection,  pertaining  to  the 
main  arc,  are  based  on  thf  input  spectrum  shown  in  Figure  2-14  which,  as 
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Figure  2-14.  Incident  electron  flux  assumed  for  auroral  arc  of  March  27,  1973 
(Black  Brant  18.205-1). 


we  have  mentioned,  is  a reduction  by  a factor  of  2 from  the  reported  values. 
However,  subsequent  comparison  of  our  results  with  data,  together  with  recent 
information  on  the  possible  importance  of  aerosol  extinction  in  the  lower 
atmosphere,  leads  us  now  to  suspect  that  the  measured  particle  data  may 
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have  been  correct  and  the  3914  photometer  data  in  error.  Although  the 
transmission  losses  due  to  aerosol  extinction  in  Alaska  at  the  time  of 
interest  to  us  are  not  known,  the  estimates  by  Photometries  (Reference  5), 
suggest  transmission  factors  at  Fort  Yukon  (due  to  Rayleigh  scattering 
and  aerosol  extinction)  of  0.60  and  0.36  for  zenith  angles  0°  and  60°, 
respectively.  This  would  then  place  the  ground-base^  optical  data  more  in 
line  with  the  full  particle  flux  data  shown  in  Figure  2-12.  As  we  will 
see  later,  other  data,  especially  those  associated  with  the  electron  density 
and  the  side- looking  photometers,  tend  also  to  support  the  higher  flux 
values.  Time  constraints,  however,  have  prevented  us  from  using  the  larger 
particle  flux  to  redo  the  calculations  reported  below.  The  reader  should 
keep  these  remarks  in  mind,  although  they  will  be  pointed  out  again  at 
appropriate  places  in  the  text  below. 

Figure  2-15  shows  the  energy  deposition  contours  for  the  auroral 
model  assumed.  The  contours  bear  a superficial  resemblance  to  the  emission 
(3914  A and  5577  X)  profiles  shown  by  Romick  and  Belon  (Reference  2Q)  for 
other  auroral  arcs,  but  they  are,  in  fact,  quite  different.  Romick  and 
Belon  subtract  out  the  "background"  in  presenting  their  data  while  we 
include  it  (the  drizzle)  in  Figure  2-15.  From  our  earlier  discussion,  it 
is  clear  that  the  magnitude,  character,  and  even  spatial  extent  of  the 
"background"  is  closely  related  to  the  main  arc,  and  a proper  interpretation 
of  the  aurora]  event  demands  that  a consideration  be  given  to  the  entire 
system  — main  arc  and  drizzle.  If  we  had  used  a more  realistic  spatial 
distribution  of  flux  in  the  drizzle  region  (rather  than  the  constant  mean 
one)  and  then  subtracted  out  a constant  background,  we  would  have  obtained 
contour  profiles  of  the  Romick  and  Belon  type. 
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Figure 


2-1 5.  Calculated  contours  of  constant  energy  deposition  inside  and 
outside  the  arc  in  magnetic  meridian  plane. 
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Figure  2-16  shows  the  calculated  energy  deposition  rate  profile 
inside  the  arc.  Shown,  for  comparison,  are  the  corresponding  results 
obtained  by  Lockheed  (LMSC) . As  mentioned  earlier,  the  differences  between 
the  two  results  can  be  attributed  almost  exclusively  to  differences  in  the 
assumed  elastic  scattering  cross  section. 

Figure  2-17  shows  the  calculated  secondary  electron  production 
rate  spectra  at  various  altitudes  in  the  arc,  while  Figure  2-18  shows  the 
secondary  electron  flux  spectra  for  the  corresponding  altitudes.  The 
secondary  flux  at  150.5-km  altitude  is  very  similar  to  that  at  122.5  km 
and  is  not  shown  in  Figure  2-18.  The  dip  that  appears  in  the  secondary 
flux  at  about  2.5  eV  is  due  to  vibrational  excitation  of  N2  that  serves 
to  deplete  electrons  with  energy  near  that  value. 

Energy  Partition  in  Arc 

The  ARCTIC  code  has  been  used  to  compute  the  volume  excitation 
rate  profiles  for  nearly  200  states  of  i\’2,  0z,  0,  CO2  and  their  ions. 
Figures  2-19  through  2-22  show  selected1  results  for  a few  of  these  spates. 

Calculational  Procedure  with  Respect  to 
Chemistry  and  Delayed  Emissions 

With  the  exception  of  the  foregoing  prouuction  rate  profiles  for 
the  main  arc,  the  results  presented  in  this  subsection  are  formatted  in  a 
way  to  facilitate  direct  comparison  with  the  rocket-  and  ground-based  data. 

To  do  this,  it  is  necessary  to  take  into  account  the  geometrical  and 
dynamical  nature  of  the  arc  as  well  as  the  rocket  position  and  motion  rela- 
tive to  the  arc.  This  means  that  because  of  arc  motion,  different  points 
along  sight  paths  (through  the  arc)  that  are  not  field  aligned  will  have  been 
irradiated  by  the  electron  flux  for  different  lengths  of  time.  Thus,  in 
computing  the  radiance  along  a given  sight  path  through  the  arc,  it  is 
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Fiqure  2-16.  Calculated  energy  deposition  rate  for  auroral  arc  of  March  27,  1973 
(Black  Brant  18.205-1)  using  incident  flux  shown  in  Figure  2-14. 


ELECTRON  ENERGY  (eV) 

Figure  2-17.  Secondary  electron  production  rate  spectra  (Black  Brant  18.205-1) 


Brant  18.205-1). 


necessary  to  allow  for  different  chemistry  times  at  each  point  along  the 

O 

path  (except  for  those  prompt  emissions,  such  as  3914  A that  do  not  depend 
on  chemistry).  For  example,  in  our  model,  when  the  main  arc  is  moving  with 
constant  velocity  to  the  north,  the  environment  at  a point  just  in  front  of 
the  arc's  leading  edge  is  that  provided  by  our  30  minute  drizzle.  On  the 
other  hand,  the  environment  at  a point  on  the  arc's  trailing  edge  corresponds 
to  that  produced  by  the  arc  flux  (+  chemistry)  operating  for  a time  AW/V  , 
where  AW  is  the  horizontal  arc  width,  and  V is  the  arc  velocity.  For 
a diffcx’ent  point,  at  a distance  Ax  behind  the  leading  edge  of  the  arc, 
the  chemistry  operates  for  a time,  t^  , given  by 


t . = Ax/V 

chem  arc 


(2-3) 


In  general,  for  an  arbitrary  point  x,y  behind  the  leading  edge,  where  x 
is  the  displacement  north  from  Fster  Dome  and  y is  the  altitude,  Ax 
is  given  by 


Ax  * x - [85  + (At  - 91)  V - (y  - 108. 7)/tan  6]  . (2-4) 

arc 

Here,  6 is  the  magnetic  dip  angle  and  At  is  the  time  elapsed  since 
rocket  launch. 

In  performing  the  chemistry  calculations  for  each  point  along  a 
given  sight  path,  care  must  be  exercised  to  see  that  the  arc  flux  is  turned 
on  while  the  point  is  within  the  arc  and  is  then  turned  off  and  the  drizzle 
flux  turned  on  when  the  arc  has  passed. 

The  results  presented  below  were  deduced  from  ARCTIC  code  runs 
that  may  be  summarized  as  follows.  The  drizzle  environment  w?^  determined 
by  inputting  the  drizzle  flux  (Figure  2-8)  at  an  altitude  of  200  km.  The 
energy  deposition  and  partition  were  then  computed  and  the  chemistry  sub- 
routine run  for  times  to  30  minutes.  The  main  arc  flux  (Figure  2-14)  was 
then  turned  on  in  the  above-determined  drizzle  environment,  and  deposition 


and  partition  profiles  were  calculated.  The  chemistry  subroutine  was  run 
for  8.2  sec  with  the  main  flux  on  and,  subsequently,  to  60  sec  with  only 
the  drizzle  flux  operating.  The  calculations,  for  chemistry  times  between 
8.2  and  60  sec,  provide  estimates  of  the  decay  in  the  emission  rates  and 
changes  in  species  concentrations  that  ensue  after  the  main  arc  has  passed 
beyond  the  point  in  question. 

Results  and  Comparisons  with  Data 

Output  from  the  ARCTIC  code,  for  both  the  main  arc  and  the  drizzle 
region,  is  voluminous.  It  includes  many  items,  such  as  species  concentra- 
tions and  volume  emission  rates,  that  will  not  be  shown  here  because  of 

* 

space  limitations.  Instead,  with  one  exception,  we  limit  the  results  shown 
to  those  that  can  be  compared  directly  with  data  available  to  us. 

Except  as  noted  below,  the  formulas  used  to  compute  the  various 
emissions  are  given  In  Reference  1 (Appendix  C)  and  will  not  be  repeated  here. 

Figure  2-23  shows  the  calculated  and  observed  zenith  radiance  at 
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3914  A as  a function  of  rocket  altitude  on  ascent.  As  was  mentioned  earlier, 
the  abrupt  change  in  slope  of  the  data  at  125  km  may  be  indicative  of  the 
presence  of  a much  larger  flux  of  electrons  with  energy  < 3 keV  than  what 
our  input  spectrum  (Figure  2-14)  shows.  As  was  also  discussed  earlier, 
the  slope  change  below  110  km  is  an  effect  related  to  the  geometry  of  the 
arc  and  to  its  motion  relative  to  the  rocket. 


The  calculated  electron  temperature  profile  is  shown  in  Figure  2-37, 
although  no  information  on  measured  values  are  yet  available 
(Reference  6) . 

The  reaction  numbering  in  the  first  column  of  Table  C-l  (Reference  1)  is 
in-  error.  From  the  top  down  it  should  read:  68,  74,  75,  71. 


57 


Figure  2-23.  3914  A brightness  viewed  vertically  upward  from 

Black  Brant  18.205-1  on  ascent. 


Considerable  weight  was  given  to  the  data  shown  in  Figure  2-23  in 
deciding  on  the  magnitude  of  the  input  particle  flux.  Use  of  the  reported 
flux  spectrum  shown  in  Figure  2-12,  rather  than  the  one  shown  in  Figure  2-14, 

o 

would  have  led  to  calculated  values  for  the  3914  A brightness  twice  as  large 
as  the  calculated  ones  shown  in  Figure  2-23  — and  more  than  twice  as  large 
as  the  indicated  data  below  130  km. 

Figure  2-24  shows  comparisons  between  the  calculated  and  observed 

O 

brightness  at  4278  A from  ground-based  scanning  photometers  at  Ester  Dome 
and  Fort  Yukon.  The  instruments  scanned  in  the  magnetic  meridian  plane, 
and  the  results,  corresponding  to  times  between  92  and  102  sec  after  rocket 
launch,  are  given  as  functions  of  the  zenith  angle. 

The  data  clearly  show  the  presence  of  a drizzle  to  the  north  and 
south  of  the  main  arc.  The  Fort  Yukon  data  give  a rough  idea  of  the  northern 
extent  of  the  drizzle  region.  The  small  peak  in  the  Fort  Yukon  data  to  the 
north  is  probably  a geometrical  effect  due  to  the  long  path  lengths  at  large 
zenith  angles. 

The  calculated  full  width  at  half  maximum  (FWHM)  of  the  arc  bright- 
ness is  narrower  than  that  measured.  We  have  performed  parametric  studies 
of  the  FWHM  seen  from  Fort  Yukon  and  Ester  Dome  for  different  arc  widths, 

AW,  and  velocities,  V . The  results  show  an  insensitivity  to  variations 
’ arc 

in  these  parameters.  We  suspect  that  the  main  source  of  broadening  is  due 
to  the  presence  of  an  "inverted  V"  structure  for  the  drizzle  region  rather 
than  the  constant  one  assumed  here. 

We  should  point  out  that  the  calculated  values  shown  in  figure 
2-24  have  not  been  corrected  for  transmission  losses  to  the  ground  if 
these  losses  are  as  great  as  quoted  earlier  (p.44),  the  magnitude  of  our 
calculated  radiances  would  be  smaller  than  those  observed.  In  fact,  if 
we  were  to  allow  for  the  variation  of  the  transmission  factor  with  zenith 
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Scanning  photometer  data  (Ref.  5)  (92-102  sec) 


angle,  and  were  to  then  increase  the  resulting  values  by  a factor  of  2,  the 
magnitude  of  our  results  at  both  lister  Dome  and  Fort  Yukon  would  appear  to 
be  in  closer  agreement  with  the  data  than  Figure  2-24  indicates.  This  would 
then  lend  support  to  the  idea  that  the  magnitude  of  the  reported  particle 

O 

flux  (Figure  2-12)  may  be  correct  and  the  rocket-borne  3914  A data  (Figure 
2-23)  erroneous. 

Figure  2-25  shows  similar  comparisons  with  the  ground-based  scan- 
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ning  photometers  operating  at  5577  A for  times  between  89  and  99  sec.  The 
calculated  volume  emission  rates  are  based  on  the  equation 

4 (5577  A)  = 1.34  [0(lS)]  (cm"3  sec"1)  (2-5) 

that  presupposes  a 0.75-sec  lifetime  for  the  JS  state  of  atomic  oxygen. 

The  O^S)  concentration  is  determined  from  the  chemistry  subroutine,  with 

the  assumed  production  mechanisms  shown  in  Table  3-2  of  the  next  section. 

* 

The  entire  chemistry  scheme  is  shown  in  Reference  1 (Table  B-l).  In  the 
next  section  we  emphasize  that  the  rate  coefficients  for  (what  now  appear 
to  be)  the  dominant  production  sources  of  O^S)  are  not  at  all  well  known. 
In  our  scheme,  the  quenching  of  N2(A3E)  by  0(3P)  is  dominant,  but  other 
reactions,  with  uncertain  rate  coefficients,  combine  to  compete  with  it. 
Direct  excitation  by  the  bombarding  electrons  appears  to  contribute  in  only 
a minor  way. 

The  foregoing  remarks  should  serve  to  indicate  that  any  agreement 

O 

between  calculated  and  observed;  radiances  at  5577  A is  probably  fortuitous. 

Figure  2-25  shows  that  the  calculated  values  are  generally  high, 
although  the  FWUM  is  approximately  correct.  The  latter  agreement  may 
imply  that  most  of  the  emission  is  from  the  main  arc.  The  angular  positions 


★ 


The  only  reaction  rate  coefficient  in  Table  3-2  that  differs  from  that  in 
Table  B-l  (Reference  19  is  for  quenching  of  N2(A3Z). 
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of  the  peaks,  however,  are  somewhat  different  from  the  measured  ones  (also 

O 

at  4278  A) . This  may  be  attributable  to  a somewhat  harder  spectrum  for 
the  main  arc  flux.  The  peak  in  the  data  to  the  north  of  Fort  Yukon  is, 
again,  probably  a geometrical  effect. 

Figure  2-26  shows  the  peak  zenith  spectral  radiance  near  5.4  pm 
observed  with  the  nose-mounted  circular  variable  filter  (CVF)  spectrometer 
on  board  the  rocket  (Reference  21)  together  with  our  computed  values.  There 
is  considerable  scatter  in  the  actual  data  points,  but  they  lie  roughly 
within  the  cross-hatched  area  shown.  We  do  not  know  what  additional  spread 
may  exist  due  to  experimental  uncertainty  associated  with  the  data  points 
themselves.  The  comparison  is  shown  for  both  the  ascent  and  descent  portions 
of  the  flight. 

The  computed  emission  is  completely  dominated  by  the  N(2D)  + 02 
chemiluminescent  reaction  for  which  photon  efficiencies  have  not  been  measured. 
However,  we  have  assumed  that  the  vibrational  states  of  NO  are  equally 
populated  up  *-0  v = 18  following  each  reaction,  and  the  photon  efficiencies 

■k 

that  result  from  this  assumption  are  shown  in  Table  C-l  (Reference  1).  On 
the  same  basis,  we  have  computed  the  spectrum  of  the  fundamental-  and  first 
overtone-band  emissions..  The  results  of  this  calculation  are  given  in 
Appendix  A. 

The  ARCTIC  code  outputs  volume  emission  rates  for  the  entire  band 
being  considered  — ir  this  case,  the  NO  fundamental  at  5.4  pm.  In  order 
to  compare  with  the  spectrometer  data,  we  used  the  information  supplied  to 
us  (References  4,  21)  for  this  instrument,  together  with  the  computed 
spectrum  shown  in  Figure  A-l.  The  wavelength  region  scanned  by  the  instru- 
ment was  from  1.6  pm  to  a maximum  of  about  5.4  pm.  \t  5.4  pm. the  bandwidth 
of  the  spectrometer  is  about  0.18  pm.  Thus.,  from  Figure  A-l,  the  peak 
emission  observed  with  this  instrument  should  occur  in  the  wavelength  in- 
terval from. 5.22  to  5.40  pm.  The  calculated  spectrum  (Figure  A-l)  shows 
that  16  percent  of  the  total  band  energy  is  emitted  in.  this  interval. 

For  correction  to  this  table,  see  footnote  on  p.  55  of  this  report. 
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Therefore,  to  compare  our  calculated  emission  rates  with  the  spectrometer 
data,  we  have  multiplied  them  by  a factor  0 . 16/ J . 18  - 0.9. 

Unless  recalibration  subsequently  changes  the  observed  values 
downward,  agreement  with  calculations  is  quite  poor,  especially  for  the 
descent  portion  of  the  flight  where  at  90  km  we  appear  to  be  about  an  order 
of  magnitude  low.  It  is  interesting  that  the  measured  radiance  on  descent 
(outside  the  arc)  is  not  much  lower  than  that  on  ascent,  even  though  a more 
pronounced  arc  enhancement  would:  be  expected  on  ascent.  Use  of  the  reported 
particle  flux  (Figure  2-12)  rather  than  one  half  of  it  (Figure  2-14)  would 
have  improved  the  agreement,  but  not  by  a sufficient  amount.  The  assumed 
photon  efficiency  for  the  N(2D)  + O2  reaction  could  be  raised,  but  not 
by  more  than  a factor  of  2 (to  preserve  energy  conservation) . It  is  more 
likely,  indeed,  that  the  value  used  here  is  already  on  the  high  side.  The 
assumed  emission  spectrum  (Figure  A-l)  may  also  put  too  much  energy  at  wave- 
lengths beyond  detectability  of  the  instrument.  However,  the  slopes  of  the 
observed  ar.d  computed  curves  are  different.  The  effect  is  as  though  an 
additional  strong  production  source  of  vibrational ly  excited  NO  (or  of 
N(2D))  is  operating  between  about  100  and  110  km  both  inside  and  outside 
the  arc.  Assuming  the  correctness  of  the  data,  we  are  presently  unable  to 
account  for  the  large  differences. 


Figure  2-27  is  another  comparison  at  5.4  pm,  this  time  with  data 
from  the  rocket-borne  dual  channel  radiometer  in  the  rocket  nose.  To  make 
the  comparison  with  the  data,  we  multiplied  our  code  values  (for  the  entire 
NO  fundamental  band)  by  a factor  f = 0.23.  This  factor  was  obtained  by 
multiplying  our  normalized  spectrum,  I (A1)  , (Figure  A-l)  by  the  instrument 
relative  response  function,  R(A),  (Reference  4)  and  integrating.  Thus, 
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Fxcept  at  altitudes  below  110  km  on  rocket  ascent,  our  computed 
values  are  again  considerably  lower  than  those  observed.  Note  that  these 
data,  like  those  from  the  spectrometer,  show  descent  radiances  nearly  as  high 
as  those  measured  on  ascent,  contrary  to  our  calculations  and  expectations. 

Figure  2-28  compares  the  calculated  values  for  the  peak  zenith 
radiance  near  4.5  pm  with  those  measured  by  the  CVF  spectrometer.  The 
theoretical  values  include  contributions  from  NO+  chemiluminescence,  CO2 
vibraluminescence,  and  ambient  (earthshine  excited)  C02  emission.  The 
N0+  and  C02  (vibraluminescence)  curves  in  Figure  2-28  were  computed  with 
the  ARCTIC  code;  the  ambient  C02  curve  was  deduced  from  the  work  of 
Kumer  and  James  (Reference  23) . 

Details  concerning  the  N0+  chemiluminescence  model  in  ARCTIC 
can  be  found  in  Reference  1 and  in  Appendix  B of  this  report.  Calculations 
based  on  this  model  show  that  most  of  the  NO+  chemiluminescence  arises 
from  the  reactions 

02+  + NO  02  + NO+ 

N2+  + 0 + N + NO+ 

that  are  sufficiently  exothermic  to  populate  the  tenth  vibrational  state  of 
N0+.  With  the  assumption  that  each  state  up  to  v = 10  is  populated  with 
equal  probability,  we  have  computed  a spectrum  for  the  emitted  photons. 

The  results  are  shown  in  Appendix  B (Figure  B-l). 

The  spectral  resolution  of  the  spectrometer  at  4.5  pm  is  about 
0.16  pm  (Reference  4).  As  the  spectrometer  sweeps  in  wavelength  over  the 
spectrum  shown  iii  Figure  B-L,  the  maximum  power  it  will  detect  will  occur 
in  a 0,16-pm  interval  from  4.28  to  4.44  pm,  repi’esenting  about  40  percent 
of  the  total  band  emission.  In  order  to  compare  our  calculations  with  the 
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Figure  2-28.  Peak  zenith  radiance  near  4.3  yin  (Black  Brant  18. 205-1 ). 
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spectrometer  observations,  therefore,  we  have  multiplied  them  by  0.40/0.16 
= 2.5,  which  then  places  our  radiances  on  a "per  micrometer"  basis. 

As  for  the  C02  vibraluminescence  (4.3-pm  band),  we  have  assumed 
a spectral  profile  similar  to  that  for  thermal  emission  at  300  °K.  Our  low 
resolution  (0.02  pm)  approximation  to  the  computed  spectrum  of  Malkmus 
(Reference  24)  is  shown  in  Figure  2-29.  As  the  spectrometer  sweeps  across 
this  spectrum,  with  a 0.16-pm  resolution,  virtually  all  of  the  band  power 
will  be  detected.  Therefore,  to  compare  our  C02  vibraluminescence  values 
with  the  spectrometer  data,  we  have  multiplied  them  by  1/0.16  =6.2. 

The  solid  curve  in  Figure  2-28  is  a sum  of  the  three  contributions 
described  above.  No  corrections  for  radiation  transport,  self  absorption, 
resonance  scattering,  etc.,  have  been  made,  although  we  don't  believe  these 
could  serve  to  enhance  the  calculated  values. 

Frcm  Figure  2-28  we  see  that  the  total  calculated  radiance  on 
ascent  is  larger  than  that  for  the  descent  portion  of  the  rocket  flight,  but 
not  nearly  as  much  as  the  data  show.  Although  the  comparisons  between  105 
and  112  km  (where  the  data  end)  on  ascent  are  quite  favorable,  the  computed 
values  below  100  km  are  low  by  as  much  as  a factor  of  4.  Considering  the 
many  uncertainties  inherent  in  the  calculation,  this  may  not  be  bad,  but 
the  rapid  fall  off  in  the  ascent  data  above  about  98  km  is  more  than  we 
can  account  for.  The  descent  data  (only  available  below  110  km):,  however, 
are  in  reasonably  good  agreement  with  the  calculations. 

Figure  2-30  is  a comparison  between  our  calculated  zenith  bright- 
ness at  4.3  pm  and  data  (Reference  21)  from  the  rocket-borne  radiometer  on 
both  ascent  and  descent.  There  is  considerable  spread1  in  the  data  that  lie 
approximately  within  the  cross-hatched  areas.  The  theoretical  values  again 
include  contributions  from  the  three  sources  described  above  and;  indicated 
in  Figure  2-30. 
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To  compare  the  code  values  (for  the  entire  bands)  with  the 
radiometer  data,  we  have  folded  the  instrument  relative  response  function 
(Reference  4)  with  the  respective  emission  spectra  (Figures  B-l  and  2-29) 
according  to  Equation  2-6.  The  resulting  correction  factors  applied  to  the 
calculated  N0+  chemiluminescent  and  C02  emissions  are  0.16  and  0.79, 
respectively.  With  these  corrections,  the  total  calculated  values  are  seen 
to  be  lower  than  the  reported  ascent  data  by  factors  ranging  from  about  6 
at  100  km  to  4 at  1:10  km  to  20  or  so  at  130  km.  The  comparison  on  descent 
is  only  slightly  better.  It  is  interesting  to  note  that  our  (full)  band 
radiances,  without  applying  the  above-mentioned  correction  factors,  would 
have  put  us  in  much  better  agreement  with  the  data. 

A peculiar  feature  of  the  radiometer  data  shown  in  Figure  2-30 
is  Its  apparent  constancy  above  about  120  kin  rather  than  a monotonic 
decline  as  would  be  expected.  The  spectrometer  data  (Figure  2-28)  are  in- 
adequate to  confirm  this  because  they  do  not  extend  to  altitudes  above  about 
112  km. 


An  increase  by  a factor  of  2 in  the  assumed  particle  flux  into 
the  arc  would  certainly  improve  our  agreement  with  the  ascent  portion  of 
this  flight,  although  it  would  be  insufficient  by  itself  to  account  in  full 
for  the  differences. 

We  turn  now  to  comparisons  with  the  side-looking  instruments  on 
board  the  rocket.  These  consisted  of  five  visible  photometers  and  a dual- 
channel radiometer  operating  in  the  5.4-  and  2.7-pm  bands.  Each  instrument 
was  mounted  with  its  viewing  axis  10°  above  the  normal  to  the  rocket  axis. 
If  the  rocket  attitude  were  vertical  at  all  times,  the  viewing  sight  paths 
would  be  10°  above  the  horizontal.  However,  as  mentioned  earlier,  the 
rocket  axis  was  inclined  about  3°  south  (magnetic)  from  the  vertical.  As 
the  rocket  spun,  therefore,  the  sight  paths  to  the  south  and  to  the  north 
were  elevated  about  7°  and  13°,  respectively,  from  the  horizontal. 
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The  data  supplied  to  us  (Reference  6)  included  uncalibrated  strip 
charts  of  voltage  output  as  a function  of  time  for  each  of  the  instruments, 
and  calibrated  versions  of  these  charts  (except  for  the  5.4-um  channel)  for 
four  specific  rocket  spins  corresponding  to  altitudes  between  95  and  103  km 
on  ascent.  Because  of  rocket  spin,  the  data  are  continuous  plots  of  the 
radiance  as  functions  of  azimuth  and  elevation  angle.  From  the  four  cali- 
brated spins  available  to  us,  we  have  attempted  to  pick  out  radiance  values 
corresponding  to  sight  paths  toward  magnetic  north  (through  the  arc)  and 
toward  magnetic  south  (away  from  the  arc)  and  have  performed  our  calcula- 
tions for  these  same  sight  paths. 

Before  comparing  the  data  and  calculations,  we  emphasize  again 
that  the  sight  paths  in  question  include  long  segments  through  the  drizzle 
region  and  (except  for  the  east-west  directions  not  considered  here)  rela- 
tively short  segments  through  the  main  arc.  In  fact,  the  radiances  measured 
to  the  south  exclude  the  arc  emissions  entirely,  at  least  until  the  rocket 
entered  the  main  arc  at  about  108-km  altitude.  Since  details  of  the  drizzle 
region  are  largely  unknown,  we  can  not  , a priori,  expect  our  assumed  model 
to  provide  agreement  in  evci’y  detail  with  the  observations.  Rather,  the 
hope  is  that  the  model  can  account,  at  least  qualitatively,  for  the  observed 
features . 


Figure  2-31  shows  our  computed  brightness  at  2.7  pm  for  the  north- 
ward- and  southward- viewing  paths.  Only  one  north-viewing  data  point  and 
two  south-viewing  ones  are  shown  because  so  many  spikes  appear  in  the  records 
it  is  hard  to  pick  out  representative  values.  In  addition,  it  should  be 
mentioned  that  the  side- looking  radiometers  (2.7  and  5.4  pm)  detected  large 
signals  that  are  believed  to  be  thermal  emission  from  ejected  pieces  of 
the  rocket  (such  as  "clamshells"  and  other  instrument  covers).  For  this 
reason  we  do  not  presently  have  data  from  these  radiometers  that  we  can, 
unequivocally,  attribute  to  the  atmosphere.  However,  the  data  points 
shown  in  Figure  2-31  are  believed  to  be  auroral ly  produced. 
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2.7-pm  brightness  (^horizontal)  viewed  toward  magnetic  north 
and  south  from  Black  Brant  18.205-1  on  ascent. 
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The  computed  emission  arises  almost  exclusively  from  the 
N(2D)  + 02  reaction.  For  comparison  purposes,  the  total  band  radiance 
(calculated)  has  been  reduced  by  a factor  of  0.63  to  allow  for  the  limited 
system  response.  This  correction  was  determined  from  Equation  2-6,  using 
the  computed  spectrum  (Figure  A-2)  and  the  relative  response  function 
reported  in  Reference  4. 

The  results  suggest  that  our  calculated  values  between  90  and 
100  km  are  low  by  more  than  an  order  of  magnitude. 

Figure  2-32  siiows  the  corresponding  results  calculated  for  the 
5. 4-pm  band.  They  have  been  corrected  for  the  instrument  relative  response 
as  described  above  (see  description  of  Figure  2-27).  Corresponding  data 
for  this  band  are  not  yet  available  to  us. 


Figure  2-33  compares  the  calculated  and  observed  radiance  at 

O 

3914  A.  In  this  figure,  as  well  as  in  the  succeeding  four  figures,  we  have 
attempted  to  extrapolate  the  four  calibrated  data  points  (for  both  the  north 
and  south-viewing  sight  paths)  to  higher  altitudes  (dashed  curves)  using  the 
uncalibrated  strip  charts  mentioned  above.  Although  this  procedure  leads, 
perhaps,  to  the  correct  trend,  it  should  be  realized  that  the  magnitudes 
shown  by  the  dashed  curves  arc  by  no  means  certain. 

In  Figure  2-33,  the  "hump"  in  the  calculated  northern  radiance 
between  95  and  110  km  is  due  to  the  main  arc.  A rather  similar  hump  appears 
in  the  corresponding  data  although  larger  in  magnitude  by  a factor  of  about 
2.  This  may,  indeed,  be  evidence  to  support  the  larger  (by  a factor  of  2) 
measured:  flux  values  (Figure  2-12)  for  the  main  arc,  instead  of  those 
(Figure  2-14)  actually  used  in  the  calculation.  Another  alternative, 
however,  is  the  presence  of  an  "inverted-V"  type  of  drizzle  that  would 
serve  to  broaden  and  increase  the  magnitude  of  the  calculated  "hump". 
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The  calculated  brightness  above  120  km  arises  primarily  from  long 
path  lengths  through  the  drizzle  region,  and  the  results  would  appear  to 
be  in  reasonable  agreement  with  the  extrapolated  data. 
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The  5577  A comparisons  shown  in  Figure  2-34  are  not  so  favorable, 
especially  above  115  km  where  the  extrapolated  data  appear  to  decrease  faster 
with  increasing  rocket  altitude.  This  may  be  evidence  for  the  presence  of 
too  many  electrons  with  energy  < 3 keV  in  our  assumed  drizzle  spectrum.  The 
"hump"  in  the  calculated  northern  brightness  is  again  due  to  the  main  arc 
and  is  in  reasonably  good  agreement  with  the  data.  As  we  have  mentioned 
before,  however,  the  production  rate  of  0('S)  atoms  is  quite  uncertain 
because  the  main  mechanisms  involved  in  the  present  calculations  are 
chemical  and  quenching  reactions  whose  rate  coefficients  are  little  more 
than  guesses. 


The  5199  A comparisons  shown  in  Figure  2-35  present  an  enigma  to  us. 
The  calculated  brightness  curves  are  based  on  the  assumption  that  the 

■k 

emission  is  from  N(2D),  with  the  volume  emission  rate  given  by 

l (5199)  = 1.23  x 10' 5 [N  (ZD) ] . (2-7) 

The  portion  of  the  sight  path  to  the  north  that  intersects  the  main  arc 
occurs  at  altitudes  sufficiently  low  so  that  the  N(2D)  should  be 
severely  quenched  by  reaction  with  02  (to  form  NO).  As  a result,  almost 
all  the  calculated  radiance  arises  from  long  paths  through  the  high-altitude 
portions  of  the  drizzle  region.  That  is  why  our  calculated  curve  (north) 
is  so  flat  and  shows  no  arc  enhancement.  By  contrast,  the  (north- viewing) 
data  clearly  show  a "hump"  at  about  100  km,  somewhat  similar  to  that 

O 

exhibited -by  the  3914  A data  in  Figure  2-33.  This  result  would  appear  to 
suggest  a strong  contribution  from  the  arc. 


Most  of  the  N(2D)  atoms  arise  from  dissociative  recombination- with  NO  . 
Contrary  to  Table  B-l  (Reference  1),  the  calculations  in  this  report  assume 
that  the  nitrogen  atoms  from  this  source  are  all  in  the  2l)  st*tc. 
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Figure  2-35.  5199.A  brightness  (-^horizontal)  viewed  toward  magnetic  north 
and  south  from  Black  Brant  18.205-1  on  ascent. 
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The  calculated  brightness  to  the  south  shows  a monotonic  increase 
with  increasing  rocket  altitude  as  the  sight  path  traverses  regions  of  in- 
creasing N(20)  concentration.  Again,  by  contrast,  the  observed  (south- 
viewing)  radiance  appears  to  exhibit  the  opposite  character,  actually  de- 
creasing monotonically  with  increasing  altitude. 

O 

Although  our  model  accounts  for  the  magnitude  of  the  5199  A 
emission  seen  to  the  north  (at  least  below  about  115  km)  it  clearly  does  not 
account  for  the  slope  of  the  observed  emission  either  to  the  nor>.h  ur  the 
south.  The  observations  seem  to  imply  the  presence  of  a strong  low-altitude 
(~  100  km)  source  of  N (2D)  both  inside  and  outside  the  main  arc.  fn  this 
connection  we  note,  as  mentioned  on  p.  65,  that  such  a source  would  also 
result  in  better  agreement  with  the  observed  zenith  radiance  at  5.4  pm. 

Further  insight  into  this  problem  may  be  gained  when  the  data  for  the  descent 
portion  of  the  flight  arc  available. 

With  respect  to  the  5199  A radiation,  we  have  also  investigated 
the  possibility  that  it  may  arise  partly  from  sources  other  than  N(2D). 

O 

The  photometer  bandpass  covers  khe  wavelength  interval  from  about  5195  A to 

o o 

5212  A (Reference  4).  Permitt  4 transitions  at  5201.8,  5199.5,  and  5206.7  A 
from  N,  N+,  and  0+,  respecti' ely , fall  in  this  wavelength  interval  (Refer- 
ence 25).  These  emissions  can,  presumably,  be  excited  by  dissociative  exci- 
tation and/or  dissociative  ionization.  The  0+  emission  could  also  be 
excited  by  direct  ionization  of  atomic  oxygen.  The  strongest  of  the  above 
three  emissions  is  likely  to  be  the  first  one  at  5201.8  A (2p25d2P  2p25p2S°) 

from  atomic  nitrogen  because  the  energy  of  the  upper  state  involved  :(d3.98  eV) 
is  considerably  lower  than  that  for  the  ot  >er  two  emissions.  In  this  case 
the  total  excitation  energy  for  the- -dissociative  excitation  process  would 
be  23.74  cV.  The  cross  section  for  dissociative  excitation  of  this  emission 
has  net  been  measured.  However,  if  we  use,  as  a guide,  the  measured  disso- 

O 

ciative  excitation  cross  section  in  Nz  for  11=77  A emission  (threshold  at 
23  eV)  we  find,  from  the  ARCTIC  code:,  brightness  curves  about  1.5  x 10  2 of 
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those  shown  in  Figure  2-33  for  the  3014  A emission.  This  would  then  provide 

O 

values  in  close  agreement  with  the  observed  5199  A radiance,  especially  to 

the  north,  and  would  account  for  the  observed  "hump"  in  the  data.  At  the 

present  time,  however,  this  possibility  is  only  a conjecture,  because  exci- 

* 

tation  cross  sections  for  the  above- cited  emissions  are  not  known. 

O 

Figure  2-36  shows  the  comparisons  at  3466  A.  The  calculations 
assume  that  the  emission  arises  from  N(2P) . Unfortunately,  the  production 
mechanisms  (mainly  chemical  reactions)  for  N(2P)  are  not  well  known  nor, 
indeed,  are  the  quenching  mechanisms.  The  production  and  destruction 
mechanisms  assumed  here  are  given  in  Reference  1. 

As  in  the  case  of  N(2D),  discussed  above,  our  calculations 

O 

show  that  the  5466  A emission  from  the  side- looking  photometer  should  be 
dominated  by  the  drizzle  region,  with  only  a negligible  contribution  (along 
the  northern  path)  coming  from  the  main  arc.  In  the  southern  direction,  our 
agreement  with  the  low-altitude  data  is  quite  good  although,  at  higher  alti- 
tudes, the  larger  calculated  values  may  imply,  igain,  an  overestimate  of  our 
assumed  drizzle  flux  for  electron  energies  < 3 kcV.  The  northern  data  seem 

o e 

to  exhibit,  the  same  type  of  "hump"  seen  at  5199  A-  and  3914  A although, 
perhaps,  not  quite  so  pronounced.  It  is  not  clear  whether  this  hump  arises 
from  an  increased  production  rate  of  N ( 2 P)  in  the  main  arc  or  from  an 
"inverted-V"  type  of  drizzle. 

o 

We  have  also  considered  the  possibility  that  the  3466  A emission 

may  receive  contributions  from  other  sources.  The  photometer  bandpass  covers 

0 ° 

the  region  from  about  5450  A to  5480  A (Reference  4).  Included  in  this  region 

+ O O O 

are  permitted  transitions  from  0 (3471  A),  Ar  (3461  A),  Ar  (3454  A),  the 

(4,6)  band  at  3457  A of  the  O2  Ilerzberg  T system,  the  (10,10)  band  at  5466 
of  the  C>2  + (I-  NFG)  system,  and  the  (3,4)  band  at  5469  A of  the  1^2(2  POS) 
system.  The  molecular  bands  appear  to  be  much  too  weak  to  play  a part, 


Furthermore,  auroral  spectrograms  do  not  clearly  show  the  presence  of  a 
line  at  5201.8  A. 
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Figure  2-36.  3466A  brightness  ( horizmtil)  viewed  toward  magnetic  north  and 

south  from  Black  Brant  18.205-1  on  ascent. 
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although  it  is  conceivable  that  the  atomic  lines  could  contribute.  Cross 
sections  for  their  excitation  are  not  available  and  we  have  not  attempted  to 
estimate  their  possible  magnitude.  However,  the  possibility  of  a signi- 
ficant contribution  from  one  or  more  of  these  lines  cannot  be  ruled  out  at 
the  present  time. 

O 

The  comparison  at  3800  A,  shown  in  Figure  2-37,  is  the  final  one 
for  the  side- looking  instruments.  The  calculated  brightness  is  based  on 
emission  from  the  (0,2)  band  of  the  N2 (2  POS)  system.  The  shapes  of  the 
observed  and  theoretical  curves  are  roughly  similar,  although  our  calculated 
values  are  low  by  factors  ranging  from  about  2 to  5.  Again,  use  of  the  full 
particle  flux  spectrum  for  the  arc  (Figure  2-12)  would  have  improved  the 
agreement. 


Figure  2-38  shows  the  calculated  electron  temperature,  T , as  a 
function  of  rocket  altitude.  The  values  at  altitudes  between  about  108  and 
140  km  correspond  to  those  in  the  arc;  other  valyes  are  thqse  appropriate 
to  the  assumed  drizzle  region.  Data  are  not  yet  available  for  comparison. 
Also  shown  in  Figure  2-38  is  the  adopted  profile  of  the  neutral  gas 
temperature. 

If  the  full  measured  arc  flux  (Figure  2-12)  had  bt^n  used  as  input 

to  the  ARCTIC  code,  the  calculated  values  for  T in  the  arc  would  have  been 
’ e 

larger  than  those  shown  in  Figure  2-38. 

A comparison  between  the  calculated-  ion  concentrations  on  rocket 
ascent  and  the  mass  spectrometer  data  of  Narcisi  (Reference  16)  is  shown  in 
Figure  2-39.  The  jagged  appearance  of  the  data  is  modulation  due  to  rocket 
spin;  their  absolute  values  were  obtained  by  normalization  to  the  election 
density  profile  deduced  from  the  plasma  frequency  probe  (Reference  16). 
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Figure  2- 37 . 3800A  brightness  (~ horizontal ) viewed  toward  magnetic  north 

and  south  from  Black  Brant  18.205-1  on  ascent. 


Figure  2-38.  Calculated  electron  temperature  at  rocket  (Black  Brant  18.205-1)  location  on  ascent. 
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Before  describing  the  comparisons  shown  in  Figure  2-39,  we  empha- 
size again  the  fact  that  the  computed  concentrations  are  based  on  chemistry 
times,  at  each  point,  determined  by  the  effective  duration  of  the  electron 
flux.  This  time,  in  turn,  depends  on  the  geometry  and  velocity  of  the  arc 
as  well  as  on  the  type,  strength,  and  prior  history  of  the  drizzle.  The 
electron  density,  in  fact,  is  rather  sensitive  to  these  parameters  which,  as 
we  have  pointed  out  before,  are  not  at  all  well  established. 

Consider  now  the  computed  and  observed  electron  density  profiles 
shown  in  Figure  2-39.  Below  about  112  km,  the  agreement  is  reasonably  good. 
Between  about  115  and  150  km,  however,  the  theoretical  values  are  lower 
than  the  measured  ones  by  as  much  as  a factor  of  2;  above  152  km  they  exceed 
the  measured  values.  As  we  noted  earlier,  energy  deposition  above  about 
120  km  occurs  mainly  by  electrons  wi jh  energies  < 3 keV.  This  is  precisely 
the  region  where  particle  flux  data  are  absent.  Our  assumed  primary 
spectrum  (Figure  2-14)  for  this  region  is  a pure  guess.  The  measured 
electron  density  profile  above  115  km  in  the  arc  (Figure  2-39)  suggests  that 
the  arc  flux  spectrum  should  increase  toward  lower  energies  below  3 keV 
(in  Figure  2-14)  rather  than  decrease  as  we  have  assumed.  On  the  other  hand, 
the  rather  abrupt  change  in  slope  of  the  measured  profile  (unlike  that  for 
the  calculated  curve)  at  about  145  km,  after  the  rocket  has  left  the  arc, 
suggests  that  we  have  overestimated  the  flux  of  low-energy  primary  electrons 
in  the  drizzle  spectrum  (Figure  2-8). 

In  summary,  several  factors  combine  to  :ender  calculation  of  the 
ion  and  electron  concentrations  uncertain,  especially  above  about  120  km. 
These  include  (1)  the  unknown  incident  flux  of  electrons  with  energy  < 3 keV, 
(2)  the  unknown  history  of  prior  electron  bombardment,  and  (3)  the  uncertain 
geometry  and  dynamics  of  the  arc  itself. 
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The  lower  (calculated)  values  of  the  electron  density  above  about 
112  km  in  the  arc  may  also  suggest  that  the  full  flux  values  (Figure  2-12) 
should  have  been  used  rather  than  the  reduced  (by  a factor  of  2)  ones 
actually  employed.  This  would  have  increased  the  calculated  values  of 
electron  density  by  about  /2,  although  the  resulting  values  near  108  km 
would  then  be  nearly  a factor  of  2 larger  than  those  measured. 

The  foregoing  remarks  are  also  apropos  to  a description  of  the  ion 
concentrations.  The  "kinks’1  that  occui  in  the  calculated  N+,  N2  + , and  0+ 
profiles  at  — 1-10  km,  as  the  rocket  leases  the  arc,  are  also  seen  in  the 
data,  but  they  then  occur  at  a somewhat  higher  altitude  and  are  not  so 
pronounced.  This  is  probably  due  to  th?  spatial  variation  of  the  drizzle 
near  the  main  arc. 

Figure  2-40  compares  the  [N0+] / [0j  + ] ratio  on  rocket  ascent. 

The  calculated  values  on  rocket  descent  in  the  drizzle  region  are  also  shown. 
The  comparison  shows  that  the  calculated  ratios  are  somewhat  higher  than  the 
measured  ascent  values. 


It  can  be  shown  (see,  for  example,  Reference  15)  that,  in  equili- 
brium, the  [NO+]/[02+]  ratio  is  dominated  by  the  quantity  ki  [N0]=/'(k?,n  ) 
where  nc  is  the  electron  density,  and  k-i  and  k?  are  rate  coefficients 
for  the  reactions 


and 


02+  -+  NO  02  + N0+ 


N0+ 


+ c *>  N + 0 


respectively.  Larger  values  for  nQ  and/or  smaller  values  for  [NO],  there- 
fore, would  serve  to  lower  our  ratios  more  in  line  with  the  data.  As  we  have 
discussed  earlier,  use  of  the  full  particle  flux  (Figure  2-12)  would,  indeed, 
increase  the  calculated  values  for  ne  and  improve  the  agreement.  As  for 
the  NO  concentration,  large  uncertainties  are  inherent  in  its  determination 
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Figure  2-40.  [NO  ]/[02]  ratio  as  a function  of  rocket  (Black  Brant  18.205-1) 
altitude;  data  and  calculations  compared. 


(as  previously  described)  because  its  magnitude  depends  on  the  prior  bombard- 
ment history  (see  Figure  2-10).  The  assumed  '“initial"  profile  of  [NO] 
prior  to  onset  of  the  drizzle  (Figure  2-11)  is  somewhat  arbitrary.  A dif- 
ferent  profile  would  have  resulted  in  different  values  for  the  [NO  ] / [O2  ] 
ratio . 

The  sensitivity  of  the  [N0+]/[02+]  ratio  to  the  strength  and  dura- 
tion of  the  bombardment  flux  is  illustrated  in  Figure  2-41.  This  figure 
shows  how  the  ratio  varies  with  time  at  an  altitude  of  103.5  km  under  the 
following  conditions:  After  50  minutes  of  continual  drizzle  bombardment, 
the  arc  flux  is  turned  on  at  t = 0 sec  and  remains  on  for  8.2  sec.  At  8.2 
sec  the  arc  is  turned  off  but  the  drizzle  remains  on.  The  results  can  be 
described  as  follows.  At  t = 0,  the  arc  flux  ionizes  02  and  drives  the 
[N0+]-/[02  + ] ratio  down  from  its  initial  value  of  about  55  to  a low  of  less 
than  2 wnen  charge  transfer  with  NO  brings  it  back  to  a stabilized  arc 
value  of  about  6.  When  the  arc  flux  is  removed,  [C_  ] decreases  by  charge 
transfer  with  NO  faster  than  the  N0+  can  decrease  by  dissociative  re- 
combination and  the  ratio  is  then  driven  to  values  exceeding  100.  Although 
not  shown  in  Figure  2-41,  the  ratio  must  eventually  (~  10  sec)  settle  back 
to  its  steady  state  drizzle  value  of  about  55. 

In  view  of  the  uncertainties  associated  with  the  bombardment 
history  of  the  region,  agreement  much  better  than  that  shown  in  Figure  2-40 
can  hardly  be  expected. 


Summary  and  Conclusions 


A summary  of  our  work  and  conclusions,  based  on  a study  of  the 
auroral  event  of  27  March  1973  near  Poker  Flat,  Alaska,  is  presented  below. 
It  should  be  emphasized  that  the  data  upon  which  much  of  this  work  was  based 
have  not  yet  been  formally  published,  but  have  been  made  available  to  us  on 
a preliminary  (subject  to  change)  basis.  Consequently,  the  conclusions 
should  likewise  be  considered  preliminary,  pending  final  confirmation  of 
the  data. 


The  summary  and  conclusions  are  as  follows: 

(1)  The  event  consisted  of  a thin  arc  (~5  km  thick)  with  IBC 
coefficient  at  least  11+ , imbedded  in  a surrounding  drizzle  region  approxi- 
mately 850  km  in  north-south  extent.  The  arc  moved  erratically  north  and 
south  in  a region  between  Ester  Dome  and  Fort  Yukon,  irradiating  the 
atmosphere  in  a manner,  ill -defined,  except  in  the  limited  air  space  sur- 
rounding the  rocket  trajectory. 

(2)  A model  to  describe  the  event,  including  the  arc  geometry, 
dynamics,  and  accompanying  drizzle  environment,  was  devised.  The  ARCTIC 
code  was  run,  for  this  model  environment,  and  comparisons  were  made  with 
the  preliminary  data. 

O 

(3)  Data  from  the  upward-viewing  rocket-borne  3914  A photometer 
are  inconsistent  with  the  particle  flux  measurements.  Either  the  particle 
flux  is  too  high  by  a factor  =2,  or  else  the  photometer  data  are  low  by  a 
corresponding  amount.  Ground-based  photometer  measurements,  uncorrected 
for  transmission  losses,  support  the  rocket  photometer  data.  However,  if 
the  ground-based  data  are  corrected  for  Rayleigh  scattering  and  aerosol 
extinction,  the  weight  of  evidence  shifts  more  in  favor  of  the  particle  flux 
data.  Comparison  with  data  from  other  rocket-borne  instruments  (especially 
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the  side- looking  3914  A and  3800  A photometers)  lends  additional  support  to 

o 

the  particle  data  and  suggests  that  the  (zenith)  3914  A rocket  data  are 
in  error. 


(4)  The  calculated  zenith  radiance  near  4.3  pm  is  about  four  times 
lower  than  the  (ascent)  CVF  spectrometer  data  below  100  km,  although  the 
agreement  is  reasonably  good  for  rocket  altitudes  between  105  and  112  km  on 
ascent  and  also  for  altitudes  down  to  90  km  on  rocket  descent.  Agreement 
with  the  (zenith)  radiance  data  from  the  dual  channel  radiometer  (4.3  pm), 
however,  is  not  so  good,  with  the  observed  brightness  above  120  km  about  20 
times  larger  than  the  calculated  values.  'Hie  apparent  constancy  of  the  radiom- 
eter data  above  120  km  with  increasing  rocket  altitude  remains  unexplained. 

(5)  The  zenith  radiance  observed  near  5.4  pm,  with  both  the  rocket- 
borne  CVF  spectrometer  and  the  dual  channel  radiometer  is  not  significantly 
different  on  ascent  (through  the  arc)  and  descent  (outside  the  arc).  This 
contrasts  with  the  calculated  results  that  show  a much  higher  relative 
brightness  from  the  arc  than  from  the  drizzle  region  outside  the  arc. 
Furthermore,  the  observed  brightness  is  about  an  order  of  magnitude  larger 
than  the  theoretical  values.  The  data  seem  to  suggest  the  presence  of  a 
source  for  vibrationally  excited  NO,  both  inside  and  outside  the  arc,  much 
greater  than  the  calculated  N(2D)  + 02  reaction  can  provide. 

(6)  Data  at  5.4  and  2.7  pm  from  the  side- looking  rocket-borne 
radiometer  are  contaminated  with  thermal  emission  from  ejected  instrument 
covers.  However,  the  detected  radiance  at  2.7  pm  that  is  (believed  to  be) 
auroral ly  produced,  is  at  least  an  order  of  magnitude  larger  than  our 
theoretical  values.  This  appears  to  be  consistent  with  the  discrepancy 
noted  in  Item  5 above. 

O 

(7)  Rocket  data  at  5199  A are  hard  to  reconcile  with  the  calcu- 
lated N(2D)  emission.  Suggested  possibilities  are  (i)  an  additional 
strong  source  of  N(2D)  operating  at  low  altitudes  (~100  km)  inside  and 
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outside  the  arc  (which  would  also  be  consistent  with  Items  5 and  6 above), 
(ii)  contributions  from  direct  excitation  of  permitted  line  radiation  (N, 

N , 0+)  that  falls  within  the  bandpass  of  the  photometer.  Dissociative 
excitation  of  such  emissions  appears  plausible. 

O 

(8)  The  shape  of  the  observed  (northward)  3466  A radiance  curve 

is  also  somewhat  difficult  to  reconcile  with  the  calculated  N(2P)  emission. 
Additional  sources  for  N(2P)  within  the  arc,  or  contributions  from  other 
sources,  are  suggested. 

(9)  Comparison  with  the  observed  electron  density  profile,  and 

O 

also  with  the  zenith  radiance  at  3914  A,  implies  that  our  assumed  particle 
spectrum  at  energies  < 3 keV  contains  too  few  electrons  in  the  main  arc  and 
too  many  electrons  in  the  drizzle  region. 

02t(1A)  EMISSION 


Enhancement  of  Oz^A)  emission  (1.27  pm)  in  auroras  has  been 
the  subject  of  controversy  for  several  years.  No  data  at  1.27  pm  are 
available  from  the  ICECAP  program.  However,  we  have  made  a comparison 
between  the  calculated  emission  from  Oz^A)  for  the  auroral  arc  of  21 
March  1973  with  observations  reported  in  the  literature  by  Gattinger  and 
Vallance  Jones  (Reference  26)  for  a different  auroral  event  that  occurred 
on  27  January  1973. 


Although  the  two  events  are  different,  one  having  occurred  near 
Poker  Flat,  Alaska,  the  other  at  Fort  Churchill,  Manitoba,  the  Poker  Flat 
event  was  evidently  more  intense  and  so,  presumably,  for  equal  deposition 
times,  it  would  give  rise  to  at  least  as  much  emission  from  O^A)  as 
would  the  weaker  event.  It  is  instructive  to  see  if  this  premise  is  indeed 
correct  and  whether  electron  bombardment  can  account  for  the  observed  en- 
hancement of  1. 27-um  emission. 
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Observations  were  made  (toward  the  zenith)  from  the  ground  at 
5577  X and  at  1.27  pm  into  IBC  Class  11+  activity.  The  5577  X emission  was 
irregular,  but  can  be  described  approximately  as  having  a peak  intensity  of 
about  50  KR  near  onset  and  a subsequent  exponential  decay  with  a time  con- 
stant ~15  minutes.  For  comparison,  we  note  that  the  observed  (and  calcu- 
lated) zenith  intensity  at  5577  X for  the  Poker  Flat  arc  was  about  70  KR. 
Therefore,  if  we  perform  our  calculations  for  the  Poker  Flat  art.  and  allow 
the  flux  to  decay  in  the  exponential  manner  described  above,  the  results 
should  overestimate  the  1.27-pm  enhancement  produced  by  precipitating 
electrons  for  the  Fort  Churchill  event. 

The  significant  processes  controlling  the.  concentration  of  02(1A) 

•k 

above  about  100  km  are  believed  to  be  electron  excitation 

e + 02(X3Z)  -*■  e + 02(a1:A)  (2-8) 

and  radiative  deexcitation 

02,(0* A)  •+  0.2  (X3Z)  + hv  . (2-9) 

Collisional  quenching  by  O^A)  with  a rate  constant  of  about  2.4  x 10  18cm3 
sec"1  should  not  be  important  above  about  90  km.  Thus,  ignoring  any  trans- 
port effects,  the  concentration  of  02.'(:1A)  is  determined  by  the  simple 
equation 

[O2 C1  A) ] + A [02(*A)]  = q(t)  . (2-10) 

Here,  A is  the  radiative  decay  rate  (2.5  * lO-4  sec  *),  and  q (t)  is  the 
volume  excitation  rate  (cm-3  sec  J)  of  02(1A). 


★ 


Chemical  excitation  has  also  been  proposed  (Reference  27). 
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In  accordance  with  observations,  we  take 

q = q0  e'at  . (2-H) 

The  solution  to  Equation  2-10  is  then 

[OzC'A)]  = (e"at  - e'At)  . (2-12) 

A decay  time  of  15  minutes  for  the  source  implies  a = 7,7  x 10  4 sec  1 . 
With  this  value  for  a,  Equation  (2-14)  then  provides  a peak  concentration 
for  02  (* A)  (and  emission  rate  at  1.27  yin)  about  35  minutes  following 
auroral  onset.  This  is  consistent  with  the  observations  of  Gattinger 
and  Vallance  Jones  (Reference  26). 

The  volume  emission  rate  in  the  (0,0)  band  at  1.27  pm  is  given  by 

4>(1.27  pm)  = A1;[02(1A)]  (2-13) 

where  Ay  is  the  transition  rate  for  the  (0,0)  band  (2.55  x 10-4  sec  *). 
Thus,  the  intensity  seen  from  the  ground  is 

OO 

1(1.27  pm)  = n J HI. 27  pm)  dz 

0 00 

= (e'At  - e~at)f  (1°  dz  (2-14) 

0 

where  n is  the  transmission  factor. 


Calculations  by  Evans  et  al.,  (Reference  28)  as  well  as  estimates 
based  on  the  calculations  of  McClatchey  and  Selby  (Reference  29)  suggest 
that  n 585  0.05.  From  the  ARCTIC  code,  we  also  find 


(2-15) 


Equation  2-14,  together  with  the  parameter  values  cited  above 
lead  to  the  results  shown  in  Fife«re  2-42.  Shown  for  comparison  are  the 
observations  of  Gattinger  and  Vallance  Jones.  The  measured  maximum  ground 
enhancement  at  1.27  pm  is  seen  to  be  about  5 times  larger  than  that  calcu- 
lated. The  two  curves  can  be  brought  into  approximate  agreement  if  it  is 
assumed  that  the  transmission  factor  is  > 0.25.  However,  the  weight  of 
evidence  seems  to  preclude  such  a large  value  for  n. 

Although  our  calculation  was  not  based  on  the  precise  environ- 
mental conditions  prevailing  at  the  time  of  the  ground  measurements,  the 
results,  nevertheless,  suggest  that  electron  bombardment  alone  is  insuffi- 
cient to  account  for  the  magnitude  of  the  observed  enhancement  at  1.27  pm. 
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Figure  2-42.  Auroral  enhancement  of  02(1a)  emission  (7.27  ym)  ; ground-based 
observations  and  calculations  compared. 


SECTION  3 

STUDIES  RELATED  TO  LWIR  PRODUCTION 
AND  TO  EUV  EFFECTS  IN  THE  AURORA 


Part  of  the  work  performed  under  the  present  contract  has  been  a 

continuation  of  the  efforts  initiated  previously  (Reference  1)  with  respect 

to  effects  produced  in  the  atmosphere  by  ultraviolet  (UV)  and  extreme 

o 

ultraviolet  (EUV,  A < 1000  A)  radiation  generated  by  auroral  electrons. 

The  motivation  for  these  studies  was  twofold.  First,  if  trapping  of  a 
sizeable  fraction  of  the  EUV  radiation  in  high- lying  states  occurs,  as  has 
been  suggested  by  Zipf  (References  30  to  33) , then  leakage  of  energy  into 
the  LWIR  might  be  significant.  In  particular,  we  sought  to  explore  the 
possibility  that  such  leakage  could  account  for  the  spectral  band  features 
observed  near  8,  9.2,  and  12  pm  in  ICECAP  and  other  auroral  experiments. 
Second,  if  the  chemistry  and  emission  properties  of  auroral  arcs  are  signi- 
ficantly altered  by  absorption  of  these  short  wavelength  radiations,  as 
has  also  been  suggested  by  Zipf,  then  a proper  interpretation  of  ICECAP 
and  other  auroral  data  would  have  to  make  allowance  for  such  effects. 

A brief  account  of  the  work  performed  and  conclusions  reached  in 
these  two  areas  is  now  presented. 

LWIR  EMISSION  FROM  ELECTRONIC  TRANSITIONS  IN  N2 

Our  efforts  to  account  for  the  LWIR  features  (near  8,  9.2,  and 
12  ym)  observed  at  high  altitudes  in  certain  rocket  experiments  (References 
34,  35)  have  centered  upon  a study  of  the  high-lying  states  of  the  N2 
moTeculo.  It  was  felt  that  N2,  being  the  dominant  molecular  species,  would 
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be  the  most  likely  source  of  CUV  and,  possibly,  of  the  observed  LWIR  bands 
that  might  arise  from  transitions  between  clcsely  adjacent  states  at  high 
energies. 


The  states  of  N2  that  we  surveyed  are  described  and  tabulated 
in  Reference  1 (Tables  5 and  6).  Over  200  states  at  energies  above  12.4  eV 
were  assembled,  mainly  from  UV  absorption  data.  Differences  were  taken 
between  those  states,  connected  by  allowed  transitions,  that  might  give  rise 
to  the  observed  LWIR  features.  A considerable  number  of  possibilities  were 
found,  although  the  uncertainty  in  wavelength  of  many  of  them  is  large  be- 
cause the  uncertainty  in  energy  of  the  levels  involved  is  often  as  great  as 
the  enei’gy  of  the  LWIR  photons  themselves. 

The  conclusions  we  have  come  to,  based  upon  a further  study  of  the 
states  involved,  are  as  follows: 

(1)  With  the  possible  exception  of  the  (0,1)  band  of  the 

w*d  >■  a 1 ITg  system  (9.5  pm)  no  states  at  energies  below  about  12  eV  appear 
capable  of  generating  the  observed  bands.  Even  the  above  band  is  usually 
weak  or  absent  under  dc  discharge  conditions  in  N2  (Reference  36)  and  is, 
therefore,  an  unlikely  candidate. 

(2)  Most  of  the  states  observed  in  UV  absorption  studies  that  lie 
above  about  12  eV,  but  below  the  first  ionization  limit  at  15.6  eV,  are 
strongly  predissociated  (Reference  37).  Those  states  that  are  observed 

in  emission  (and  hence  not  predissociated)  do  not  have  suitably  adjacent 
states  with  which  to  connect  to  produce  the  observed  LWIR  bands. 

(3)  For  states  lying  above  the  first  ionization  potential,  high 

resolution  absorption  studies  have  shown  that  the  rotational  lines  involved, 

0 

although  varying  in  width,  are  of  the  order  of  0.01  to  0.1  A wide  (Doppler 
width  ~0.002  A).  This  corresponds  to  lifetimes  of  ~10-11  to  10~12  sec  and 
implies  that  the  states  are  strongly  autoionized  and/or  predissociated 
(Reference  38) . 
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If  leakage  of  trapped  EUV  energy  into  the  LWIR  is  to  account  for 
the  observed  features,  then  not  only  should  suitably  adjacent  states  exist 
in  the  molecular  system,  but  also  the  cross  sections  for  UV  absorption  into 
and  radiation  from  them  should  be  large.  The  foregoing  observations, 
however,  appear  to  preclude  the  satisfaction  of  these  conditions,  and  we 
conclude  that  the  N2  molecule  is  not  the  source  of  the  LWIR  features. 

The  attempts  outlined  above  to  explain  the  LWIR  bands  were  predi- 
cated on  the  assumption  that  they  are  aurorally  produced.  However,  observa- 
tion of  the  features  under  aurorally  quiet  conditions  above  Poker  Flat, 

Alaska  in  February  1974  suggests  that  this  may  not  be  the  case.  Consequently, 
no  further  study  of  the  bands  was  made  pending  clarification  of  the  circum- 
stances surrounding  their  occurrence. 


UV  GENERATION  AND  ITS  EFFECTS  ON  5577  A EMISSION 

The  generation  of  UV  and  EUV  radiation  by  precipitating  auroral 
electrons,  and  the  subsequent  effects  upon  the  chemical  and  optical  pro- 
perties of  auroras,  is  not  well  known,  to  say  the  least.  Good  laboratory 
and  field  data  on  the  EUV  spectrum  are  much  needed  for  this  purpose. 
References  39  to  41  give  some  information,  but  it  is  far  from  adequate. 

The  expected  data  from  the  University  of  Pittsburgh  are  not  yet  available. 
Consequently,  firm  conclusions  regarding  the  importance  of  the  effects 
cannot  yet  be  made. 

We  have,  nevertheless,  performed  a cursory  investigation  of  certain 
aspects  of  the  problem.  First,  an  attempt  was  made  to  estimate  the  fraction 
of  precipitated  electron  energy  that  may  go  into  EUV  radiation  by  studying 
the  energy  partition  calculated  using  the  ARCTIC  code.  This  is,  aair.i t+edly;, 
rather  crude  because  the  cross  sections  for  exci-tation  of  the  higher- lying 
states  of  the  major  species  are  not  well  known  nor,  indeed,  is  the  ultimate 
fate  of  the  excitation  energy.  Second,  we  have  considered  the  influence  of 
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OUV  radiation  on  the  excitation  of  0(!S)  and  the  subsequent  emission  at 

O 

5577  A.  Tli is  green- line  emission  from  atomic  oxygen  is  a standard  reference 
in  auroral  work,  including  the  ICECAP  program,  and  it  is  important,  therefore, 
to  understand  the  mechanisms  involved  in  its  generation. 

A brief  summary  of  the  results  of  this  work  are  as  follows: 

Estimated  Fraction  of  precipitated  Electron  Energy 
Radiated  at  A < 1000  A. 

In  order  to  estimate  the  fraction  of  precipitated  electron  energy 
radiated  in  the  EUV  above  100  km,  the  ARCTIC  code  was  used  to  calculate  the 
energy  partition  resulting  from  the  stopping  of  primary  electrons  in  an 
air  mixture  consisting  of  73%  N2 , 14%  02,  and  13%  0.  This  mixture  corres- 
ponds to  an  altitude  of  about  120  km. 

The  partition  into  187  different  states  was  determined  for  primary 
electron  energies  ranging  from  1 eV  to  100  kcV.  The  detailed  results  are 
not  presented  here,  although  a somewhat  condensed  version  for  an  N2/02 
mixture  is  shown  in  Table  2-1. 

The  history  of  many  of  the  states  following  their  excitation  is 
not  well  known.  States  of  N2  and  02  that  predissociate  and  preionize 
may  ultimately  lead  to  radiation  from  the  product  atoms  or  atomic  ions. 

The  fraction  of  such  energy  radiated  and  its  spectrum  are  presently  unknown. 
Fortunately,  however,  only  a relatively  small  part  of  the  incident  electron 
energy  is  deposited  in  states  such  as  these. 

Table  3-1  is  our  best  guess  of  the  fraction  of  energy  radiated  in 
the  EUV  and  the  approximate  wavelengths  involved.  It  is  based  on  the 
partition  determined  for  the  stopping  of  a 10  keV  electron,  although  the 
detailed  partition  is  actually  rather  insensitive  to  the  primary  energy  for 
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values  above  about  1 keV.  Most  of  the  radiation  is  seen  to  arise  from 
Rydberg  transitions  in  N2  and  from  the  b'1!!*  state  of  N2.  The  indicated 
fraction  radiated  by  the  latter  state  is  probably  an  upper  limit  because  some 
of  the  higher  vibrational  states  associated  with  it  are  not  observed  in 
emission  and  are  evidently  predissociated.  Direct  excitation  of  atomic 
oxygen  appears  to  contribute  only  negligibly  to  the  EUV. 


The  conclusion  from  the  foregoing  (rather  crude)  analysis  is  that 
not  more  than  about  10  percent  of  the  precipitated  electron  energy  is 
radiated  in  the  EUV. 


Table  3-1.  Estimated  fraction  of  incident  electron 
energy  (10  keV=)  radiated  in  the  EUV. 


Wavelength  (A)  Radiating  Species 


Process 


Fraction  Radiated 


900-970 

N2(b  '£*) 

direct  excitation 

1 .7" 2 

~ 960 

N2( Rydberg) 

direct  excitation 

6.0"2 

~ 865 

N2 (Rydberg) 

direct  excitation 

1.8'2 

916 

N+ 

dissoc.  ioniz.  of  N2 

1.4"3 

776 

N+ 

dissoc.  ioniz.  of  N2 

7 . 7“** 

747 

N+ 

dissoc.  ioniz.  of  N2 

5.6"4 

672 

N+ 

dissoc.  ioniz.  of  N2 

6.3"“ 

533 

N+ 

dissoc.  ioniz.  of  N2 

4.9"1* 

879 

0 

dissoc.  excit.  of  02 

9. O'5 

990 

0 

dissoc.  excit.  of  02 

1.8~4 

833 

0+ 

dissoc.  ioniz.  of  02 

2.1-4 

717 

0+ 

dissoc.  ioniz.  of  02 

9.2"5 

TOTAL 

^l.O'1 

1 imit) 
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Possible  Effects  of  EUV  radiation  on  5577  A Emission 

Two  possible  effects  of  EUV  radiation  on  the  5577  X emission  from 
atomic  oxygen  were  considered.  One  deals  with  the  peculiar  observations 
reported  by  Romick  and  Belon  (Reference  3-13)  that  shows  that  the  thickness 

O 

of  an  auroral  arc  measured  at  5577  A is  only  about  half  as  large  as  that 
measured  at  3914  X.  The  other  deals  with  the  importance  of  O^S)  excita- 
tion by  photodissociation  of  O2  as  compared  to  other  excitation  mechanisms. 

Arc  Width  at  5577  X 

Auroral  excitation  of  O^S)  and  subsequent  emission  of  the 

O 

5577  A line  is  normally  considered  to  arise  from  a combination  of  direct 
plus  chemical  processes.  The  ARCTIC  code,  for  example,  predicts  that  nearly 
half  of  the  excitation  arises  from  the  quenching  reaction 

N2(A3Z*)  + 0(3P)  + N2(XV)  + 0(*'S).  (3-1) 

Most  of  these  processes,  however,  including  Reaction  3-1,  should  produce  an 
excitation  rate  nearly  proportional  to  the  energy  deposition  rate.  This 
means  that  the  emission  rate  from  0(1S)  should  be  nearly  proportional 
(except  for  lifetime  considerations)  to  that  from  N2  at  3914  X which  arises 
from  direct  electron  excitation.  Therefore,  photometers  viewing  upward 
along  the  magnetic  field  lines:  and  scanning  across  an  auroral  arc,  should 
see  about  the  same  intensity  fall -off  from  the  center  to  the  edge  of  tie 
arc  at  both  5577  A and  3914  X.  The  results  of  Romick  and  Belon,  hoV’cver, 

O 

show  that  the  brightness  at  5577  A falls  off  from  the  arc  center  faster 
than  does  that  at  3914  X. 

According  to  Zipf  (Reference  32) , the  dominant  excitation 
mechanism  for  0(1S)  in  the  aurora  is  photodisspciation  of  02  by 
aurorally-generated  EUV.  We  therefore  sought  to  explore  the  possibility 
that  trapping  of  EUV  radiation  near  the  arc  center,  with  subsequent  photo- 
dissociation  of  02,  might  account  for  the  RomicI:  and  Belon  result. 
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A detailed  treatment  of  radiation  trapping  in  an  arc  is  clearly 
complex  and  involves  a consideration  of  radiation  transport.  We  attempted, 
instead,  to  determine  the  plausibility  of  the  idea  by  estimating  the  (rela- 
tive) energy  density  of  EUV  radiation  across  an  arc  based  on  a one-dimensional 
random  walk  model  for  the  photons.  In  particular,  we  assumed  a gaussian 
profile  for  the  volume  production  rate  of  EUV  photons  across  the  arc.  Each 
photon  was  then  allowed  to  random  walk  in  one  dimension  with  a probability 
ai  of  resonant  scatter  and  (1  - a i ) of  absorption  at  each  scattering.  Thus, 
if  the  volume  production  rate  of  EUV  photons,  as  a function  of  distance  x 
from  the  arc  center,  is  given  by 

P(x)  = P(0)  e“x2/2x°2  , (3-2) 


it  can  then  be  shown  that,  for  the  simple  model  considered,  the  photon 
number  density  across  the  arc  is  determined  from  the  equation. 


n.  = 


_ P(0) 


nccL 


ai 


m-i 


mAo  ‘ 


Xo 


x0 


exp 


2 (mA0 2 + xoz) 


)• 


(3-3) 


m=i 


Here,  n is  the  molecular  concentration,  c the  speed  of  light,  aT  the 
total  scattering  cross  section,  and  Ao  the  radiation  mean  free  path. 


Equation  3-3  is  a superposition  of  gaussian  functions  each  of 
which  is  broader  than  the  volume  production  rate  profile  3-2.  Consequently, 
the  energy  density  of  the  radiation  field  falls  off  more  slowly  from  the 
arc  center  than  does  the  photon  production  rate.  But  the  photon  (EUV) 

O 

production  rate  should  be  proportional  to  the  production  rate  of  3914  A 
photons  since  both  are  proportional  to  the  energy  deposition  rate.  But  if 

O 

the  5577  A emission  is  excited  mainly  by  (EUV)  dissociative  excitation 
of  02,  then  its  volume  emission  profile  should  behave  like  that  of  the  EUV 
energy  density.  The  result  should  then  be  a broadening  of  the  arc  width 

O O 

seen  at  5577  A relative  to  that  seen  at  3914  A,  contrary  to  the  observations 
of  Romick  and  Be Ion. 
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The  conclusion,  subject  to  the  limitations  of  the  simple  model 
adopted,  is  that  the  Romick  and  Belon  observations  cannot  be  accounted  for 
through  photodissociative  excitation  of  O^S)  by  trapped  EUV  photons. 

Photo-Dissociative  Excitation  of  0(*S) 

The  quantum  yield  of  O^S)  atoms  resulting  from  the  process 

02  + hv  0 + 0(rS)  (3-4) 

has  been  measured  by  Lawrence  and  McEwan  (Reference  42) . For  wavelengths 

O O 

outside  the  range  850  A < A < 1200  A the  yield  was  found  to  be  negligible. 

O 

The  peak  yield  occurs  for  A ~ 1050  A and  attains  a value  of  about  10  percent. 

O 

This  means  that  about  10  percent  of  the  photons  (A  ~ 1050  A)  that  are 
absorbed  by  02  are  effective  in  producing  0(JS)  atoms.  The  average 
value  over  the  indicated  wavelength  range  is  about  5 percent. 

An  upper  limit  to  the  volume  excitation  rate  of  O^S)  atoms  by 
process  3-4  can  now  be  estimated  as  follows.  Let 

<f  = energy  deposition  rate  by  bombarding  electrons 

o o 

Hi  = fraction  of  energy  radiated  in  interval  850  A to  1200  A 
r\z  = fraction  of  the  radiated  energy  absorbed  by  02 
ri3  = quantum  yield  of  0(JS)  atoms 
The  excitation  rate  is  then  given  by 

jd[0dt"S)Jjuv  * ^2113  <*?/  hv  (3-5) 

where  hv  is  the  average  photon  (UV)  energy. 
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Reliable  estimates  for  m are  not  available.  However,  study  of 

our  results  on  energy  partition  in  an  N2,  O2,  0 mixture  (mentioned  above) 

* 

leads  us  to  assign  a value  of  about  0.14.  This  is  very  likely  an  upper 
limit.  Other  upper  limits  are  02  = 1>  and  = 0.1.  For  an  average  photon 
wavelength  of  1000  A (12.4  eV) , Equation  3-5  then  leads  to  the  result 

jlMSHj  < , x 10-=  j . (3-6) 

I dt  luv 

As  an  example,  consider  the  auroral  arc  of  27  March  1973  described 
in  Section  2.  At  the  altitude  of  peak  energy  deposition  («  101  km),  we  have 
£ = 2.7  x 10 7 eV  cm"3  sec"1.  Equation  3-6  then  yields 

jd[Q^S)}j  < 2 ? x 1q4  (cm-3  sec“i)  . (3-7) 

( dt  )UV 

Because  of  the  upper  limits  assumed  for  the  parameters,  especially 
rj 2 » it  seems  likely  that  the  value  given  by  Equation  3-7  is  too  large  by  at 
least  a factor  of  2.  Nevertheless,  it  is  instruct/ive  to  compare  it  with 
values  for  the  dominant  excitation  processes  included  in  the  ARCTIC  code. 
Table  3-2  shows  the  calculated  excitation  rates  for  0(1S)  from  these  other 
processes  at  an  altitude  of  101  km  in  the  auroral  arc,  8.2  sec  after  onset. 
With  the  assumed  rate  coefficients  (none  of  which  is  well  known)  nearly 
half  of  the  excitation  rate  comes  from  quenching  of  N2  (A3£) . Zipf  (Refer- 
ence 32)  has  suggested  a rate  coefficient  for  this  quenching  reaction  of 

1.5  x 10" 11  cm3  sec"1.  Use  of  this  value  would  provide  an  even  larger  con- 

■*  * 

tribution  from  A- state  quenching. 


* 00 

This  includes  a contribution  of  0V04  from  wavelengths  1000  A to  1200  A o 
along  with  the  0.10  contribution  from  Table  3-1  for  wavelengths  < 1000  A. 

A value  of  1.5  x 10" 11  for  this  rate  coefficient,  along  with  the  other  ex- 
citation, processes  shown  in  Table  3-2,  leads  to  5577  n intensities  larger 
than  observed  on  ICECAP.  Consequently,  the  rate  coefficient  has  been 
lowered  to  the  value  shown  in  Table  3-2. 
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Table  3-2.  Volume  excitation  rates  of  OCS)  at  101  km 
in  auroral  arc  of  2 7 March  1973. 


PROCESS 

ASSUMED  RATE 
CONSTANT 
(cm3  sec"1) 

VOLUME  EXCITATION 
RATE  _ 

(cm-3  sec"1) 

direct  electron  excitation 

6.43 

0+(2D)+02 

► o(rs)+ot 

r10 

1.83 

0+(2p)+0? 

-o(xs)+ot 

I"10 

3.62 

N+(3P)+02 

-0(1S)+N0+ 

T 1 o 

1.3“ 

N+(‘S)+02 

- 0(:-S)+N0+ 

1- 1 0 

l.l3 

n(2p)+o2— 

-0(1S)+N0 

1 - 1 1 

2.1“ 

e+02 — -^0(3P)+0.97  0(‘D) 

2.8-7 

5.43 

+0.03  0(,S) 

N2(A3E)+0(3P)- 

N2  ( X 1 Z ) +0  ( 1 S ) 

5,0" 12 

4.2“ 

£*  9-1" 


If  the  upper- limit  contribution  due  to  the  UV  (Equation  3-7)  is  now 
added  to  the  results  shown  in  Table  3-2,  we  find  that  it  would  contribute 
about  20  percent  to  the  total  excitation  rate  of  O^S).  Realistically, 
the  contribution  is  probably  no  greater  than  about  10  percent  if  the  con- 
tributions from  the  other  processes  are  reasonably  accurate.  This  result 
is  at  variance  with  the  above-mentioned  conclusions  of  Zipf  who  attributes 
most  of  the  excitation  of  OfVS)  to  the  UV  mechanism. 

109 


'1  - 


It  should  be  cautioned,  however,  that  the  only  excitation  rates 
shown  in  Table  3-2  that  we  have  some  confidence  in  are  those  arising  from 
the  direct  electron  process  and  from  dissociative  recombination  in  O2. 
Rate  coefficients  for  the  other  processes  are  very  uncertain.  Therefore, 
no  firm  conclusions  on  the  l'elative  importance  of  the  UV  mechanism  for 

o 

excitation  of  5577  A emission  can  presently  be  made,  and  our  preliminary 
conclusion  is  highly  tentative,  to  say  the  least. 


SECTION  4 


ATMOSPHERIC  HEAVE  UNDER  AURORAL  ELECTRON 
BOMBARDMENT  CONDITIONS 


The  MICE  code  has  been  used  to  estimate  the  atmospheric  heave 
resulting  from  IBC  Class  III  auroral  electron  bombardment  conditions.  The 
primary  electron  flux  at  the  top  of  the  atmosphere  was  taken  to  be 

4>  = 2.2  x 10s  e . (cm  2 sec-1  eV  J)  (4-1) 

This  corresponds  to  an  integrated  energy  deposition  rate  of  1.6  * 1013  eV 
cm  2 sec  1 . 


The  altitude  dependence  of  the  deposition  rate  (ergs  gm-1  sec-1) 
determined  by  the  ARCTIC  code  is  shown  as  the  solid  curve  in  Figure  4-1. 

For  input  to  the  MICE  code  we  have  approximated  this  curve  by  the  dashed 
curve  in  Figure  4-1,  namely,  a deposition  rate  of  1„4  x 105  (ergs  gm-1 
sec-1)  above  100  km. 

The  pressure  efficiency  above  100  km  is  not  well  known.  However, 
in  orde~  to  obtain  an  upper  limit  on  the  atmospheric  heave,  we  have  arbi- 
tral .y  assumed  an  efficiency  of  70  percent  so  that  a constant  value  of 
1 x 105  (ergs  gm-1  sec-1)  is  input  as  a heating  rate  above  100  km. 

The  results  are  shown  in  Figures  4-2  and  4-3.  Figure  4-2  shows 
the  altitude  of  selected  air  parcels  as  functions  of  time  for  the  first  4 
minutes  of  the  bombardment.  Figure  4-3  shows  the  vertical  displacement  of 
the  air  at  4 minutes  as  a function  of  the  initial  altitude.  From  Figure 
4-3  wc  see  that  the  maximum  heave  occurs  at  an  altitude  of  about  145  km 
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Figure  4-2.  Atmospheric  heave  resulting  from  auroral  electron  bombardment. 


and  attains  a vertical  displacement  of  less  than  4 km  after  4 minutes  of 
continual  bombardment.  This  is  only  about  one  fourth  of  the  scale  height 
at  that  altitude.  The  actual  amount  of  heave  is  expected  to  be  less  than 
the  results  shown  would  indicate  because  we  have  used  a pressure  efficiency 
that  is  undoubtedly  too  high. 

We  conclude  that  atmospheric  heave,  resulting  from  ele  -tron 
precipitation  during  an  auroral  event,  is  not  a significant  effec..  This 
conclusion  is  borne  out  by  recent  incoherent  scatter  measurements  during 
the  recovery  phase  of  a magnetic  substorm  (Reference  43) . The  same 
measurements  suggest,  however,  that  heave  at  high  altitudes  resulting 
from  joule  heating  by  the  electric  field  that  drives  the  auroral  electrojet, 
may  serve  to  enhance  the  F- region  ionization  by  significant  amounts. 

It  should  be  noted,  nevertheless,  that  evidence  cited  in  Section  2 
shows  that  the  auroral  zone  latitudes  undergo  almost  continual  electron 
precipitation  ranging  from  a weak  "drizzle"  to  the  intense  bombardment 
associated  with  strong  arcs.  The  calculations  described  above  pertain  to 
heave  effects  expected  from  strong  events  that  are  typically  of  short 
duration.  We  have  not  addressed  the  question  of  what  may  happen  to  the 
atmosphere  at  high  latitudes  when  it  undergoes  a "drizzle"  - type  bombard- 
ment for  days  on  end.  It  may,  conceivably,  produce  a significant  departure 
of  the  composition,  or  even  the  density,  from  that  prevailing  at  lower 
latitudes. 
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SECTION  5 

SUMMARY,  DISCUSSION,  CONCLUSIONS  AND  RECOMMENDATIONS 


The  work  reported  above  has  been  divided  into  three  main 
categories:  (1)  Applications  of  the  ARCTIC  code,  (2)  EUV  effects  in  the 
aurora,  and  (3)  Atmospheric  heave  from  auroral  electron  precipitation.  The 
largest  effort,  by  far,  has  gone  into  category  (1),  with  particular  emphasis 
on  a study  of  the  auroral  event  of  27  March  1973  near  Poker  Flat,  Alaska. 
This  section  gives  a brief  summary  of  the  results  and  conclusions  of  our 
studies  in  these  three  areas,  and  lists  some  specific  recommendations  for 
future  action. 

APPLICATIONS  OF  THE  ARCTIC  CODE 

For  potential  use  in  weapon  effects  applications,  the  ARCTIC  code 
was  first  used  to  determine  the  energy  partition  that  results  from  stopping 
of  keV  electrons  in  normal  air.  The  calculated  partition  into  the  various 
excited  states  was  tabulated.  The  accuracy  of  the  results  is  not  known 
because  many  of  the  cross  sections  employed  are  quite  uncertain.  However, 
the  computed  energy  required  to  produce  an  ion  pair  in  air  agrees  well  with 
the  limited  experimental  data. 

The  code  was  then  run  for  the  case  of  a hypothetical  auroral 
event,  and  the  resulting  energy  deposition  profile  was  compared  with  ones 
predicted  by  other  codes.  Our  results  were  found  to  be  in  good  agreement 
with  those  from  the  Lockheed  code,  that  performs  the  calculation  in  quite 
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a different  way.  Minor  differences  in  the  results  can  be  attributed  to 
differences  in  the  assumed  elastic  scattering  cross  sections.  Species 
excitation  rates  were  also  calculated  and  the  results  presented. 

The  major  application  of  the  ARCTIC  code  was  to  the  auroral  event 
of  27  March  1973,  cited  above.  With  the  aid  of  data  supplied  to  us  by  AFCRL, 
Photometries,  and  Utah  State  University,  we  first  constructed  an  approxi- 
mate model  to  describe  the  geometrical  and  dynamical  aspects  of  the  event 
that  consisted  of  a large  drizzle  region  surrounding  a smaller,  but  more 
intense,  main  arc.  ARCTIC  runs  were  made  to  assist  in  constructing  the 
model,  and  subsequent  runs  were  used  to  provide  radiance  and  other  data 
for  comparison  with  observations.  Conclusions,  based  on  these  comparisons, 
are  as  follows: 


(1)  At  the  time  of  rocket  entry  to  the  arc,  the  arc  had  a width 
a 5 km  (north-south)  and  was  moving  northward  with  a velocity  «0.61  km 
sec-1.  The  associated  drizzle  region  was  probably  of  the  "inverted-V" 
type  and  extended  for  a distance  of  about  540  km- to  the  north  and  310  km 
to  the  south  of  the  main  arc. 

O 

(2)  Rocket  photometer  data  of  zenith  brightness  (3914  A)  are 
inconsistent  with  the  particle  flux  measurement.1;.  Either  the  particle  data 
are  too  high  by  a factor  ^ 2 or  else  the  photometer  data  are  too  low  by  a 
corresponding  amount.  Although  our  calculations  were  performed  with  the 
assumption  that  the  particle  data  are  too  high,  subsequent  comparisons  and 
considerations  (including  allowance  for  atmospheric  transmission  losses 

inherent  in  the  ground-based  photometer  data)  suggest  that  it  is  probably 

o 

the  3914  A photometer  data  that  are  in  error. 

(3)  The  calculated  zenith  radiance  near  4.3  ym,  including  con- 
tributions from  N0+  and  C02 , cannot  account  for  the  reported  magnitude 
of  the  observed  emission.  Although  the  data  froim  the  spectrometer  and 
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radiometer  are  not  entirely  consistent,  our  theoretical  values  are  generally 
lower  than  those  observed.  The  apparent  constancy  of  the  radiometer  data 
with  increasing  rocket  altitude  above  120  km  (if  real)  is  a mystery. 

(4)  The  calculated  zenith  radiance  near  5.4  pm  (NO  chemilumin- 
escence) is  as  much  as  an  order  of  magnitude  lower  than  the  obseived  values 
from  both  the  dual  channel  radiometer  and  the  CVF  spectrometer.  Ti.e  meas- 
ured zenith  radiance  on  rocket  descent  (viewing  outside  the  arc)  is  nearly 
as  great  as  that  measured  on  ascent  (viewing  through  the  arc) . This  sur- 
prising result  is  not  borne  out  by  the  calculations. 

(5)  Some  of  the  data  from  the  rocket-borne  side-looking  radiom- 
eters (2.7  and  5.4  pm)  are  contaminated  by  thermal  emission  from  ejected 
parts  of  the  rocket.  However,  other  pieces  of  the  data,  believed  to  be 
auroral ly  produced,  suggest  that  our  computed  NO  radiances  (at  2.7  pm) 
are  low  by  at  least  an  order  of  magnitude.  This  result  is  consistent  with 
the  findings  in  Item  4 above. 

(6)  The  shapes  of  the  observed  (north-  and  south-viewing) 
radiance  curves  at  5199  A are  hard  to  reconcile  with  the  calculated  N(  D) 
emission.  Suggested  possibilities  include  (i)  an  additional  strong  source 
of  N(2D)  operating  at  low  altitudes  (~  100  km)  inside  and  outside  the  main 
arc  (which,  incidentally,  would  be  consistent  with  Items  4 and  5 above), 
and  (ii)  contributions  from  direct  excitation  of  permitted  line  radiation 
(N,..  N+,  0+)  that  falls  within  the  bandpass  of  the  photometer. 

O 

(7)  The  shape  of  the  observed  (northward)  3466  A radiance  curve 
is  also  somewhat  difficult  to  reconcile  with  the  calculated  N(2P)  emission. 
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(8)  Comparisons  with  the  observed  electron  density  profile,  and 

O 

also  with  the  zenith  radiance  at  3914  A,  suggest  that  our  assumed  particle 
spectrum  at  energies  < 3 keV  contains  too  few  electrons  in  the  arc  and  too 
many  electrons  in  the  drizzle  region. 

(9)  The  positive  ion  mass  spectrometer  data  imply  that  the  NO 
concentra- ion  between  100-  and  105-km  altitude  was  > 109  cm"3  in  the  vicinity 
of  the  auroral  arc.  This  represents  nearly  a 40-fold  enhancement  over  cor- 
responding low-latitude  values. 

(10)  Comparisons  between  the  calculated  zenith  radiance  at  1.27  ym, 
and  observations  reported  in  the  literature,  suggest  that  auroral  enhancement 
of  O2  C1^)  f mission  cannot  be  accounted  for  entirely  on  the  basis  of  electron 
precipitation. 

With  respect  to  the  foregoing  calculations,  we  should  mention, 
perhaps,  that  no  account  was  taken  of  the  electro jet.  The  work  of  Cole 
(Reference  44)  implies  that  at  altitudes  below  150  km,  the  effect  of  the 
electrojet  may  be  described  (by  means  of  a conductivity  tensor  in  the 
manner  of  Rees  and  Walker  (Reference  45))  as  joule  heating.  The  ion  tem- 
perature can,  conceivably,  be  raised  above  the  neutral  gas  temperature, 
although  the  unknown  character  of  the  electron  precipitation  over  the 
drizzle  region  would  render  difficult  any  reasonable  assessment  of  the 
electric  field  and  the  resulting  increase  in  T^  . However,  temperature 
increases  by  more  than  factors  of  2 or  3 would  not  be  expected,  and  we 
believe  that  our  neglect  of  the  electrojet  does  not  introduce  serious  errors. 

In  concluding  our  remarks,  we  wish  to  emphasize  that  the  calcu- 
lations and  comparisons  presented  should  be  regarded  as  a "first-cut" 
effort.  Information  obtained  from  additional  and/or  revised  data,  along 
with  insight  gained  from  the  above  comparisons,  will  place  us  in  a position 
to  define  a better  model  for  the  auroral  event  and  to  reach  more  definitive 
conclusions  in  the  foliow-on  effort. 
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EUV  EFFECTS  IN  AURORA 


Two  areas  addressed  with  respect  to  possible  extreme  ultraviolet 

(EUV)  effects  in  the  aurora  were:  (a)  the  possibility  that  aurorally- 

generated  EUV  radiation,  trapped  in  high-lying  states  of  N2 , might  leak 

out  into  the  LWIR  through  transitions  to  neighboring  states,  thereby 

accounting  for  certain  LWIR  features  observed  from  ICECAP  rockets  at  high 

o 

altitudes;  (b)  modification  of  the  5577  A emission  through  EUV  trapping 
and  subsequent  production  of  O^S)  by  photodissociative  excitation  of 
02. 


The  conclusions  reached,  which  are  very  tentative  because  of  the 
almost  total  lack  of  data  on  the  auroral ly- generated  EUV  spectrum,  are  as 
follows: 

(1)  The  N2  molecule  does  not  appear  to  be  the  source  of  the 
observed  long  wavelength  infrared  (LWIR)  features.  Bands  that  arise  from 
states  at  energies  below  12  eV  in  N2  occur  at  the  wrong  wavelengths. 
States  above  about  12  eV,  whose  energy  differences  make  them  attractive 
candidates,  are  strongly  predissociated  and/or  preionized  and  so  will  emit 
weakly,  if  at  all. 

(2)  Based  on  an  ARCTIC  code  study  of  energy  partition,  we  esti- 
mate that  not  more  than  about  10  percent  of  the  precipitated  electron-- 
energy  is  radiated  in  the  EUV. 

(5)  A simple  model  of  photon  transport  leads  to  the  conclusion 

O 

that  trapped  EUV  photons,  that  produce  5577  A emission  by  photodissociative 
excitation  of  02 , cannot  account  for  the  observed  fact  that  auroral  arcs 

o o 

are  narrower  when  measured  at  5577  A than  when  measured  at  3914  A. 
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(4)  No  more  than  about  20  percent  of  the  5577  A auroral  emission 
can  be  attributed  to  photodissociative  excitation  of  0(aS)  by  auroral ly- 
gererated  UV  radiation.  This  result,  that  is  at  variance  with  the  findings 
of  Zipf,  is  based  partly  on  the  validity  of  uncertain  rate  coefficients  for 
chemical  excitation,  incluaing  the  quenching  reaction 

N2(A3E)  + Of'P)  N2 (X1^)  + O^S), 
which  we  find  to  be  the  dominant  source  of  O^S)  atoms. 

AURORALLY-PRODUCED  ATMOSPHERIC  HEAVE 

To  study  the  possibility  of  atmospheric  heave  under  auroral 
bombardment,  the  ARCTIC  code  was  first  run  to  establish  an  altitude  profile 
of  energy  deposition  under  Class  III  auroral  conditions.  An  approximation 
to  the  calculated  profile  was  then  used  as  input  to  the  MICE  code  and, 
taking  what  we  believe  to  be  an  upper  limit  for  the  pressure  efficiency, 
we  calculated  the  heave  resulting  from  4 minutes  of  continual  bombardment. 

The  conclusion  is  that  heave,  under  the  stated  conditions,  is 
insignificant.  The  maximum  effect  occurs  at  145  km  where  the  vertical 
displacement  of  the  air  is  only  one  fourth  of  the  scale  height  at  that 
altitude. 


In  spite  of  the  foregoing  conclusion,  it  seems  possible,  however, 
that  the  effects  of  prolonged  bombardment  in  the  auroral  oval,  that  may 
occur  almost  continually,  could  lead  to  significant  differences  between 
the  high-  and  low-latitude  profiles.  Resolution  of  this  question  would 
entail  a calculation  that  takes  into  account  additional  effects  such  as 
winds,  diffusion,  solar  radiation,  etc. 


RECOMMENDATIONS 


In  order  to  improve  our  understanding  of  the  mechanisms  involved 
in  a disturbed  environment,  and  to  maximize  the  usefulness  of  future  field 
data,  the  following  recommendations  are  made.  It  is  realized  that  some  of 
them  may  be  difficult  to  implement.  They  are  listed  in  approximate  order 
of  importance  as  we  currently  view  them. 

(1)  Rocket-borne  optical  and  infrared  instruments  should  be 
oriented  so  that  they  view  upward  along  the  magnetic  field  lines.  Since 
the  particle  flux  data  permit  a determination  of  the  energy  deposition  and 
partition  along  the  field  lines,  direct  comparison  between  the  computed  and 
observed  emissions  would  then  be  facilitated. 

Instruments  that  look  sideways  receive  radiation  from  long  paths 
through  (drizzle)  regions  that  are  not  well  defined,  if  at  all,  by  the 
rocket-borne  particle  detectors.  Calculated  radiances  for  such  paths  are 
very  uncertain  and  comparisons  with  the  measured  values  are,  therefore, 
much  less  valuable. 

We  realize  that  space  limitations  probably  prevent  the  mounting 
of  all  optical/infrared  instruments  in  the  rocket  nose.  However,  if  it 
were  possible  to  orient  the  rocket  so  that  its  axis  were  nearly  horizontal, 
then  the  side- looking  instruments  could  view  upward  along  the  field  lines 
once  per  rocket  spin.* 

(2)  Rocket  Blights,  coordinated  with  satellite  measurements  of 
particle  flux  over  a wide  region  surrounding  the  arc,  would  be  very 
desirable.  This  is  especially  important  if  Item  1 above  cannot  be  imple- 
mented, because  information  is  needed  to  define  the  drizzle  region  that 
contributes  in  a dominant  way  to  the  radiance  along  sight  paths  that  are 
not  field  aligned  through  the  arc.  In  fact,  such  satellite-acquired  data 

C 

* Lacking  this  possibility,  ft  is  important  that  at  least  the  5199  A photom 
eter  be  mounted  so  as  to  view  in  the  prescribed  manner. 
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over  a period  of  days  prior  to  a proposed  event  would  help  to  define  the 
environmental  conditions  prevailing  at  the  time  of  rocket  launch. 

(3)  Instrument  covers  should  not  be  ejected  from  the  rocket,  if 
possible.  Their  presence  in  the  field  of  view  of  the  (IR)  instruments 
seriously  degrades  the  usefulness  of  the  data. 

(4)  A better  understanding  of  the  important  excitation  mechanisms 

. O 

for  0 (lS)  atoms,  leading  to  5577  A emission,  can  be  achieved  if  rate 
coefficients  for  the  following  reactions  are  measured: 

N2(A3Z)  + 0(3P)  -►  N^E)  + O^S) 

N (2P)  + 02  -*■  NO  + O^S) 

N+(3P)  + 02  -*>  N0+  + 0(>S) 

0+(2D)  + 02  -*-=02+  + 0 ( 1 S) 

These  are  the  dominant  production  sources  for  5577  A emission  currently  in 
the  ARCTIC  code,  but  their  rate  coefficients  are  quite  uncertain.  Direct 
electron  excitation  of  0(*S)  is  not  competitive  with  these  sources  nor, 
do  we  think,  is  photodissociation  of  02  the  dominant  mechanism.  However, 
these  conclusions  could  change  if  our  assumed  rate  coefficients  are  too 
high. 


(5)  Rate  coefficients  involving  the  production  and  destruction 
of  N(2P)  atoms  are  also  needed.  These  include 

nVd)-  + 02  N(2P)  + 02  + 

N+(lS)  + 02  -*•  N(2P)  + 02  + 

N(2P)  +.  02  NO  + 0 
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(6)  Photon  spectra  and  yields  from  the  following  chemiluminescent 
reactions  need  to  be  measured: 

N(2D)  + 02  ->  NO  + 0 
02  + + NO  N0+  + 02 
N2+  + 0(3P)  N0+  + N 

For  the  third  reaction,  the  branching  ratio  for  production  of  N(‘*S)  and 
N(2D)  atoms  should  also  be  obtained. 

(7)  In  situ  measurements  of  the  EUV  spectrum  in  an  auroral 
environment  are  needed.  Without  such  data,  the  importance  of  EUV  effects 
on  auroral  interpretations  will  remain  inconclusive. 


(5.4,  2.7  pm) 
(4.3  pm) 
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APPENDIX  A 


SPECTRAL  DISTRIBUTION  OF  HO  CHEMILUMINESCENT  EMISSION  FROM 
THE  REACTION  N{ 2 D)  + 02~  NO  + 0 


The  spectral  distribution  of  energy  emitted  in  the  fundamental 
(5.4  pm)  and  first  overtone  (2.7  pm)  bands  of  NO  from  the  N(2D)  + O2  re-  ■ 
action  has  not  been  measured.  It  can  be  calculated  if  the  distribution  among 
vibrational  states  of  NO  immediately  following  the  reaction  is  known.  Unfor- 
tunately, this  information  is  not  available.  However,  in  conformity  with 
the  assumptions  made  earlier  with  respect  to  chemiluminescence,  we  will 
assume  that  the  vibrational  levels  are  populated  with  equal  probability  up 
to  the  maximum  allowed  by  energy  conservation.  On  this  basis,  we  proceed 
to  derive  the  spectral  distributions,  making  due  allowance  for  the  rotational 
structure  of  the  bands. 


If  n is  the  highest  vibrational  state  permitted  by  energy  conser- 
vation (n=18  for  the  N(2D)+02  reaction),  we  then  assume  that  the  probability 
of  populating  any  of  the  levels  v=0,T,2,---,  n is  the  same.  With  this 
assumption,  it  can  then  be  shown  that,  for  the  fundamental  bands,  the  average 
number  of  transitions  (v+v-1)  that  occur  following  cascade  is  given  by 


A 


V,V-1 


= N(v) 
n+1 


\,;v-l/  (Av,v-1+\,v-2-* 


(A-l) 


where  Ay  y 1 and  \ v 2 are  radiative  transition  probabilities  for 
transitions  v-^v-1  and  v+v-2,  respectively.  The  quantity  N(v)/(n+l),  which 
is  the  average  occupation  of  vibrational  level  v per  molecule,  is  determined 
from  the  equation 


Preceding  paes  Hank 
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N(v)  = 1+N(v+1)  • V1,.v/CVl,vtVl,v-l) 
*N(V2)  • Avt2jV/Avt2>v.Avt2>v^) 


(A-2) 


with  N(n)  = 1;  N(n.l)  = N(n»2)  = 0. 


The  transition  probabilities  (sec-1)  corresponding  to  the  funda- 
mental and  first  overtone  vibrations  are,  respectively. 


A 


, = 2 . 8x10" 8 S0  v2  .•  v 
, V- 1 u v,v-l 


and 


A = 1 .4xio_ 8 S0  v2  • v(v-l) 
v,v-2  u v,v-2  v ’ 


(A-3) 

(A-4) 


1 - *>  - 1 

where  S0  and  So  are  the  integrated  absorption  coefficients  (cm  atm  ) 

at  N.T.P.  for  the  fundamental  and  first  overtone  bands,  respectively.  The 

wavenumbers  (cm-1)  of  the  fundamental  and  first  overtone  transitions  are 

given,  respectively,  by 


and 


v _ = 0)  -2v  0)  x + (3v2+1/4)  a)  y 
x v,v-2  e e e v e'e 

v = 2[w  - (2v-l)  w x + (3v2-3v+7/4)  u)  y ] 
v,v-2  1 e *•  1 e e "■  7 e7eJ 

Similarly,  for  the  first  overtone  bands, 

'/tv,v-2  = "n+1^"  Av,v-2,/ ('Av,v-l+Av,;v-2^  ’ 


(A-5) 

(A-6) 


(A-7) 


Table  A-l  shows  the  values  adopted  here  for  the  various  molecular 
constants  appearing  above  and  below. 
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Table  A-1.  Values  adopted  for  molecular  constants  of  N0(X2n). 


CO 

(cm-1) 

w * 

e e 

(cm-1 ) 

w v 
e-'e 

(cm-1 ) 

So 

(cm"2atnf 1 ) 

1 

So 

(cm_2atm"1) 

Be 

( cm' 1 ) 

ae 

(cm'1) 

A 

1904.03 

13.97 

-0.0012 

128 

2.11 

1.7046 

0.0178 

1 

Let  ffv  ,)  bo  the  fraction  of  the  total  fundamental  band 
1 v,v-l' 

energy  emitted  per  reaction  in  the  v-*v-l  transition.  It  follows  from  the 
foregoing  that 


*v,v-l  V,-V-l^^. Vv,v-1  'y*v,v-l . 

* V=1 


Similarly,  for  the  overtone  bands. 


(A-8) 


f^v,v-2)  Vv,v-2'/V,v-2/£2Vv,v-2;yv,v-2* 


(A-9) 


Table  A-2  shows  the  calculated  values  of  f for  each  of  the 
bands  comprising  the  fundamental  and  first  overtone  band  systems. 

Each  vibrational  band  within  the  fundamental  and  first  overtone 
band  systems  has  a rotational  structure  composed  of  many  rotation-vibration 
lines.  The  lines  form  part  of  the  P,  Q,  or  R branches  of  the  band  if  the 
rotational  quantum  number  changes  by  -1,  0,  or  +1,  respectively.  In  partic- 
ular, the  frequencies  (cm"1)  of  the  vibration-rotation  lines  of  the  P,  Q,  and 
R branches  of  a band  that  arises  from  a v'+v"  transition  are  given  by* 
(Reference  46) : 

*The  treatment*presented  here  ignores  fine  structure  effects  due  to  a non-zero  electron 
spin  of  the  ground  electronic  state  of-NQ.  For-the  resolution  considered,  this  is 
entirely  adequate. 
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Table  A-2.  Calculated  values  for  f(v  ) and  f(v 

V J V* * V y V~ 


Transition 

A Cum) 

f 

Transition 

A(ym) 

T-  0 

5.330 

1.47"1 

2-  0 

2.685 

2-  1 

5.411 

1.29"1 

3-  1 

2.726 

3-  2 

5.494 

1 . 1 3" 1 

4-  2 

2.768 

4-  3 

5.579 

9.86-2 

5-  3 

2.812 

5-  4 

5.667 

8.6T2 

6-  4 

2.857 

6-  5 

5.758 

7 . 51 “2 

7-  5 

2.903 

7-  6 

5.853 

6.52"2 

8-  6 

2.951 

00 

1 

5.950 

5.63"2 

9-  7 

3.001 

9-  8 

6.050 

4.83-2 

10-  8 

3.052 

10-  9 

6-.il  54 

4.12"2 

I T-  9 

3.105 

11-10 

6.261 

3. 47" 2 

12-10 

3.160 

12-11 

6.372 

2.88~2 

13-11 

3.217 

13-12 

6.488 

2 . 35“ 2 

14-12 

3.276 

14-13 

6.607 

1 .88" 2 

15-13 

3.337 

15-14 

6.731 

1.44- 2 

16-14 

3.401 

16-15 

6.859 

O 

t 

ro 

17-15 

3.467 

17-16 

6.993 

6. 81" 3 

18-16 

3.536 

18-17 

7.131 

3 . 99~ 3 

J 


5P  * Vv"  * (BV'-BV"3  j2  + A2  - 2BV. 


V'  V' 


V"  V' 


(J=0,l,2, ) (A- 10) 


vQ  = Vv"  + (®vt“Bv»5  j2  + J + ^"“V3  A*  C^0’1^,— ) 


C A- 11) 


= \tv»  + J + (B,.'B„..)  j2  + (B,.rB„)  A*  (<7*1*2, 3, — ) 


(A- 12) 

Here,  v , ,,  is  the  frequency  in  the  absence  of  any  rotational  considerations 

(given  by  Equations  A- 5 and  A- 6).,  J is  the  rotational  quantum  number  of  the 

upper  state  involved  in  the  transition.  A is  the  quantum  number  associated 

with  the  electron  angular  momentum  about  the  internuclear  axis,  and  is 

determined  in  terms  of  the  rotational  constants  B and  a by  the  equation 

6 0 


Bv  = Be  - <*e(v+l/2).  (A- 13) 

Values  for  the  rotational  and  other  constants  used  above  are  shown 
in  Table  A-l. 


If  a Boltzmann  distribution  among  the  rotational  states  is  assumed, 
the  intensities  of  the  rotational  lines  within  the  P,  Q,  and  R branches 
are  then  given  by  the  following  expressions  (Reference  46): 


K? 

Ip(V’.J)  = “Q 


issmMl 1:A>_  a'hcF'',CJykT(Jso,i,2,.-) 


n:  .i run  e ^ j 


VV',J)  = Q J(J+l) 


(A- 14) 
(A-l 5) 
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IR(v’,J)  = 


K VR  (2J ' +1)  A2 

Qr  J(J+1) 


-hcF  , (J)/kT 

e v (J=l,2, 3, — ) 


(A- 16) 


Here,  K is  a constant,  Qy  is  the  rotational  part  of  the  partition  function, 
and  Fv (J) , the  rotational  energy  (cm-1)  of  a state  with  vibrational  and  rota- 
tional quantum  numbers  v and  J,  respectively,  is  given  by 


Fy(J)  = By  J(J+1)  - BvA2.  (A- 17) 

The  foregoing  relations,  together  with  the  results  shown  in  Table 
A-2,  are  sufficient  to  determine  the  spectral  distribution  of  energy  emitted 
by  NO  from  the  N(2D)+02  reaction.  In  particular,  the  fraction,  gp(v',J),  of 
the  energy  of  a given  band  (specified  upper  (v')  and  lower  (v")  vibrational 
quantum  numbersj  that  is  emitted  in  a P-branch  transition  from  an  upper 
state  with  rotational  quantum  number  J,  is  just 


gp(V,J)  = rp(v',J)/I(v»),  (A- 18) 

where 

I(v')  = £[Ip(v*,J)  + lQ(v',J)  + IR(v’,J)].  (A- 19) 

J 

Thus,  for  example,  the  fractional  contribution  of  this  transition  to  the  total 
energy  emitted  in  the  fundamental  band  system  of  NO  is  f(\>v,  v,  j)  gp(v',J). 

In  computing  the  spectral  distributions,  we  have  included  contribu- 
tions for  rotational  states  up  to  J=  30.  The  fractional  contributions  of 
all  lines  falling  in  a given  wavelength  interval  were  then  summed.  The 
wavelength  intervals  (spectral  resolution)  chosen  for  the  fundamental  and 
overtone  band  systems  are  0.033  ym  and  0.025  ym,  respectively.  The  tempera- 
ture was  taken  to  be  250  °K,  corresponding  to  conditions  at  about  110-km 
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altitude.  The  results  are  shown  in  Figures  A-]  and  A-2.  In  these  figures, 
the  ordinate  gives  the  fraction  of  the  total  energy  emitted  in  the  band 
system  within  the  incremental  wavelength  interval  chosen.  For  example, 
from  Figure  A-l  we  see  that  5.1  percent  of  the  energy  emitted  per  reaction 
in  the  fundamental  band  system  occurs  between  5.4  and  5.433  ym. 
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Figure  A-2.  Calculated  spectrum  of  NO  (first  overtone)  emission  from  the  N(2D)+02  reaction. 


APPENDIX  B 

SPECTRAL  DISTRIBUTION  OF  NO+  CHEMILUMINESCENT 
, EMISSION  FROM  THE  REACTIONS  , 

02  + NO  02  + NO  AND  N2+  + 0 -*•  NO  + N(4S) 


The  ARCTIC  code  calculates  the  chemiluminescent  emission  from 
the  4. 3-pm  fundamental  band  of  NO+  by  including  contributions  from  the 


reactions : 

02+  + NO  -*  0?  (X32)  + N0+  + 2.8  eV  (B-l) 
N2+  + 0(3P)  * N(4S)  + N0+  + 3.10  eV  (B-2) 
N2+  + 0(3P)  ■*  N(2D)  + NO+  + 0.73  eV  (B-3) 
0+(4S)  + N2 (v=0, 1 ,7)  N(4S)  + NO+  (B-4) 


Since  the  photon  efficiency  for  each  reaction  is  presently  unknown,  the 
usual  assumption  is  made  that  one  half  of  the  available  energy  goes  into 
vibrational  excitation  of  NO+.  On  this  basis,  the  calculations  described 
in  this  report  have  shown  that  most  of  the  NO+  chemiluminescence  arises 
from  Reactions  B-l  and  B-2.  The  exothermicity  of  these  two  reactions 
is  sufficient  to  populate  the  vibrational  levels  of  N0+  up  to  v = 10. 

The  procedure  used  to  determine  the  spectral  distribution 
of  the  emitted  photons  is  similar  to  that  described  in  Appendix  A and  will 
not  be  repeated  here.  The  main  difference  between  the  treatment  here  and 
that  for  NO  in  Appendix  A arises  because  N0+  possesses  a 2 ground 
electronic  state  and,  consequently,  the  Q-  branch  lines  are  missing  from 
its  spectrum. 
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The  molecular  parameters  for  NO  are  not  well  established.  The 
values  adopted  here  are  shown  in  Table  B-l  together  with  the  corresponding 
references.  Values  for  So,  and  So',  the  integrated  absorption  intensities 
of  the  fundamental  and  first  overtone  bands,  respectively,  are  substantially 
lower  than  the  values  reported  by  Stair  et  a!.,  (Reference  49) . However, 
only  the  ratio  between  these  quantities  is  involved  in  our  calculations, 
and  for  this  ratio  the  results  reported  in  References  48  and  49  are 
reasonably  close. 

Table  B-l.  Values  adopted  for  molecular  constants  of  N0+(X! E) 


Parameter 

Value 

Reference 

to  (cm-1) 
e 

2377 

47 

Ve(cni_1) 

16 

47 

Be(cm_1) 

2.002 

47 

ae(cm_1) 

0.020 

47 

S0(cm~ 2 atm-1 ) 

88.8 

48 

S0'(cnf2  atm-1) 

0.6 

48 

The  frequencies  of  the  vibration- rotation  lines  of  the  bands  are 
determined  from  Equations  A-5,  A-6,  A-iO,  and  A-12  with  A = 0. 


The  fraction  of  the  total  energy  of  the  fundamental  band  system 
emitted  per  reaction  in  the  v -*■  v - 1 transitions  and,  correspondingly,  for 
the  first  overtone  system  in  the  v -*■  v - 2 transitions,  is  shown  in  Table 
B-2  (compare  Table  A-2  for  the  case  of  NO). 


Table  B-2.  Relative  importance  of  contributing  bands  to 
the  fundamental  and  first  overtone  systems 
in  N0+  chemiluminescence. 


Transition 

A(pm) 

f 

Transition 

X(pm) 

f 

1-0 

4.264 

2 . 02” 1 

2-0 

2.147 

6 . 23” 2 

2-1 

4.323 

1.75"1 

3-1 

2.177 

1.07'1 

3-2 

4.384 

1.50"1 

4-2 

2.208 

1.35'1 

4-3 

4.446 

1 . 26~ 1 

5-3 

2.239 

1 . 49” 1 

5-4 

4.511 

1 . 04“ 1 

6-4 

2.272 

1.50'1 

6-5 

4.577 

8.40"2 

7-5 

2.305 

1.39"1 

7-6 

4.645 

6.50"2 

8-6 

2.340 

1 .1 8” 1 

1 

00 

4.715 

4.74'2 

9-7 

2.375 

CM 

1 

LO 

00 

00 

9-8 

10-9 

4.787 

4.861 

3.10'2 
1.56' 2 

10-8 

2.412 

5.12"2 

The  spectral  distributions  were  calculated  assuming  a Boltzmann 
distribution  of  rotational  states  at  a temperature  of  250  °K.  The  wave- 
length interval  (spectral  resolution)  chosen  for  the  fundamental  and  first 
overtone  band  systems  is  0.02  pm,  The  results  are  shown  in  Figures  B-l 
and  B-2.  In  these  figures,  the  ordinate  gives  the  fraction  of  the  total 
energy  emitted  in  the  band  system  within  the  incremental  wavelength  interval 
chosen.  For  example,  from  Figure  B-l  we  see  that  about  3.1  percent  of  the 
energy  emitted  per  reaction  in  the  fundamental  band  system  occurs  between 
4.2  and1  4.22  Pm. 


Figure  B-2.  Calculated  spectrum  o 


APPENDIX  C 

LISTING  OF  THE  ARCTIC  CODE 


This  appendix  provides  a listing  (on  microfiche)  and  a brief 
description  of  the  ARCTIC  code. 

The  ARCTIC  code  has  the  schematic  structure  indicated  in  Figure 
C-l.  The  (0,0)  overlay  consists  of  the  main  program  ARCTIC  which  is  essen- 
tially a driver  program  and  three  subroutines,  ECRD,  LATS,  and  ATMS.  These 
subroutines  are  common  to  many  of  the  subsequent  overlays  called.  During 
execution,  the  ARCTIC  main  .program  calls  one  of  the  seven  overlays  CHEM, 

ELDEP,  UVDEP,  IFAC,  EDIT,  APOT,  PIPUT  depending  on  the  function  to  be 
performed.  These  functions  are  described  below. 

Overlay  (7,0)  (main  program  PIPUT)  reads  in  all  data  and  initializes 
all  common  blocks  and  large  core  memory,.  This  overlay  also  calculates  all 
storage  parameters,  initializes  the  main  deposition  grid  (if  no  restart  is 
involved)  and  determines  whether  it  is  a one,  two  or  three  dimensional 
problem. 


Overlay  (2,0)  (main  program  ELDEP)  deposits  the  energy  from  the 
incoming  electrons  and  partitions  it  among  the  different  states.  In  the 
case  of  a two-  or  three-  dimensional  problem  the  entire  deposition  is  done 
before  returning  to  the  main  program  ARCTIC. 

Overlay  (3,0)  (main  program  UVDEP)  is  at  present  a dummy  one. 

It  is  included  in  order  to  make  future  provision  fo-  the  transport  and 
deposition  of  UV  radiation  generated  by  the  bombarding  electrons. 

Preceding  page  blank 
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Figure  C-l.  Schematic  structure  of  the  ARCTIC  code. 
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Overlay  (5,0)  (main  program  EDIT)  is  an  output  routine,  which  j 

prints  all  relevant  information  on  the  electron  and  UV  depositions. 

Overlay  (4,0)  (main  program  IFAC)  is  an  interface  routine  used  : 

for  converting  data  between  the  overlays  (2,0)  (ELDEP) , (3,0)  (UVDEP)  and 
overlay  (1,0)  (CHOM) . This  routine  allows  some  flexibility  in  altering 
the  deposition  and  chemistry  overlays  separately. 

I 

Overlay  (6,0)  (main  program  APOT)  is  a dummy.  1 

I 

Overlay  (1,0)  (main  program  CHEM)  uses  the  output  of  overlay 
(4,0)  to  perform  the  subsequent  chemistry.  This  overlay  actually  performs 

the  integration  of  the  coupled,  time-dependent  differential  equations.  j \ 

\ 

. \ 

. 4 

The  microfiche  listing  of  the  code  is  enclosed  in  the  envelope  j 5 

attached  to  the  inside  of  the  back  cover.  ! ^ 

! 1 

i l 

! 1 
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